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Abstract 

The article reviews recent developments on magnetic properties of superconductors with 
anisotropic Cooper pairing. In particular, we show how the concept of broken symmetries is 
applied to the investigation of the mixed state in superconductors with a multicomponent order 
parameter. Starting from the phenomenological description in the framework of the generalized 
Ginzburg-Landau theory, we discuss different types of quantized vortices appearing at Hd in 
states with and without time-reversal breaking. General classification of superconducting phase 
transitions in a uniform magnetic field at Hc2 is constructed. Vortex lattices of different forms 
are found in the vicinity of the upper critical field. Symmetry arguments are used to classify 
phase transitions inside the mixed state. Special attention is given to results which can be 
obtained analytically. Also special emphasis is put on the open questions of the theory. 



1 Introduction 



The BCS microscopic theory of superconductivity [^] was preceded by the phenomenological ap- 
proach to the superconducting phase transition developed in the works by Ginzburg, Landau [^°] 
and Abrikosov [^]. Considering the superconducting order parameter as a scalar complex wave func- 
tion which breaks gauge symmetry, one can explain many remarkable properties of superconductors: 
the Meisner effect, critical behavior of bulk samples and thin films, two types of superconductiv- 
ity, and so on. The microscopic approach developed later clarified the physical meaning of the 
complex order parameter and established the relation between phenomenological constants of the 
Ginzburg-Landau (GL) functional and microscopic characteristics of superconducting metal. 

One of the basic points of the BCS theory is the assumption that electrons are paired in a fully 
isotropic s-wave state with both total spin S and orbital momentum L of the Cooper pairs equal to 
zero. The possibility of the unconventional superconductivity with anisotropic, non s-wave pairing 
has been discussed since early 60th. In particular, it was shown [^'^'^'^] that anisotropic Cooper pairs 
are formed if the interaction between fermion quasiparticles is attractive for at least one value of 
the relative orbital momentum. 

The feature of the phenomenological theory of unconventional superconductors is that they 
are described by the multicomponent order parameter with a more complicated structure of the 
GL functional. This leads to qualitatively new properties of unconventional superconductivity. The 
behavior of superconductors with a multicomponent order parameter in magnetic field is the subject 
of this review. 

Until the beginning of the 80th, supcrfiuid "^He was the only known example of an unconventional 
pairing. However, because of the absence of electric charge, the interaction of superfiuid component 
in '^He with magnetic field is restricted to the weak spin paramagnetism of Cooper pairs, and lies 
therefore beyond our attention. (The complete review on superfluity of ^He is given in [^^].) 

Since the discovery of superconductivity in so called heavy fermion compounds '''23,51 they 
have been extensively discussed as new possible candidates for anisotropic pairing. Some of them, 
such as CeCu2Si2, UBcia, URu2Si2, UPta, UPd2Al3, UNi2Al3, reveal various unusual properties, 
which could be attributed to unconventional superconducting states p6,66,68j 

The most convincing evidence of the nontrivial superconductivity in heavy fermions is the compli- 
cated H-T phase diagram of superconducting UPta (Fig. 1). In numerous experiments the splitting 
of the superconducting phase transition in zero magnetic field, the discontinuity in slope (kink) of 
Hc2 (T) , two different superconducting states at if = and three different vortex phases have been 
observed (see [36>68j g^^^j reference there in). Such multiphase diagram could not be interpreted in 
terms of the usual GL theory. It seems natural to explain it assuming that superconductivity in 
UPta is described by a multicomponent order parameter with different minimums of corresponding 
GL functional, and phase transitions take place between them. 

A number of questions about magnetic properties of unconventional superconductors are inspired 
by these experiments: what are the common and different features of the unconventional and s- 
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wave superconductors, what are the phase transitions in the mixed state of superconductors with 
a multicomponent order parameter? In this paper we report on the recent progress in this field 
achieved since the pubhcation of previous reviews [22.65j concentrate mostly on theoretical 

aspects of the description of the vortex state in unconventional superconductors (the list of the 
discussed questions is given at the end of Sec. 1.3). We do not suppose to present here details of the 
concrete models designed to explain peculiar properties of UPt3, this was done recently in a number 
of works [2''''^^'^^] . However, we use these models as an illustration of the general approach and try 
to stress, when it is possible, which particular properties of the mixed state of multicomponent 
superconductors can be used for reliable identification of an unconventional superconductivity in 
heavy fermions. 

Three following Subsections have an introductorily character and give the outline of the GL 
theory of unconventional superconductors, the main points of the Abrikosov theory of the mixed 
state of conventional s-wave superconductors and the properties of multicomponent superconductors 
at = 0. 

1.1 Ginzburg-Landau Theory 

In this Subsection we consider GL theory of unconventional superconductivity developed in [8,72,75j^ 
We define the order parameter and the symmetry group of the space uniform state of unconventional 
superconductors (at H = 0) which will be used to describe the mixed state in the next Sections. 
The two-component model, which is commonly used for the description of unconventional super- 
conductivity in UPt3, is considered as a particular example of unconventional superconductivity. 

The superconducting second order phase transition in the ordinary case results in the breaking 
of normal state gauge symmetry. To give the phenomenological description of the phase transition 
in an unconventional superconductor, we first define the symmetry group of a normal metal, which 
can be broken when the condensate of Gooper pairs appears. This group includes: 

• The gauge group U{1) of multiplication of electron wave functions ^ by an arbitrary phase 
factor: * e'"*. 

• The space group of the crystal lattice. 

• The group of transformations of the spin space. 

• The time reversal operation R. 
Two remarks are relevant at this point. 

(i) The space group, generally speaking, includes different combinations of translations, rotations 
and refiections. However, because of the large size of Gooper pairs (which at T = is of the order 
of coherence length ^o) in comparison to lattice constant a (in BGS theory a/^o ~ Tc/Ep), the 
discrete crystal structure has no effect on the classification of superconducting states. Therefore, 
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a superconductor can be considered as uniform but anisotropic medium with the space symmetry 
T X G, where T is the abehan group of continuous three dimensional translations, and G is the 
point group of the crystal. Classification of superconducting states were done taking into account 
discrete crystal structure in p^'^'^]. 

(ii) Transformations of the spin space are separate from real space transformations only if the 
spin-orbital interaction is negligible (as in the case of superfluid ^He). In the opposite case of 
strong spin-orbital coupling, one should assume that point group G acts in both coordinate and 
spin spaces. Since the strong spin-orbital coupling seems to be relevant for properties of heavy 
fermion compounds, we will concentrate on the latter case. Note, however, that the spin-orbital 
interaction does not lift the Kramers degeneracy, and superconducting electrons are characterized 
by the two-component wave function \E'q,. This pseudospin degeneracy is important for the triplet- 
singlet classification of the superconducting state ["^■''^jsj ^^^^ below). 

Finally, the symmetry group of the normal state Q can be written as follows: 



It is more convenient to rewrite in Wigner variables: Cooper pair center of mass coordinate 
r = (ri + r2)/2 and position on the Fermi surface k, which is obtained by Fourier transformation: 
(ri — —>■ k. The electron energy spectrum is reconstructed in a thin layer near the Fermi surface. 
Hence, the dependence of AQ,^(k, r) on the direction of k is more important than on its modulus. 
In this Subsection we assume that the superconducting state is homogeneous in space what means 
that Aq,^ does not depend on r. For this reason we will not discuss here translational properties of 
the order parameter. 

The possible superconducting states are classified according to the symmetry of Aa^(k). We 
have already mentioned that the superconducting phase transition breaks the gauge invariance U{1). 
Hence possible residual symmetries of Aa^(k) may be found by listing all subgroups of group Q which 
do not contain the group U{1) as a separate subgroup f'"^^]. For a conventional superconductor this 
subgroup is equal to G x R. If it is smaller than G x R, the order parameter Aa^(k) is defined to 
be unconventional. 

Not all subgroups of Q correspond to energeticaly stable superconducting phases. The applica- 
tion of the Landau symmetry approach restricts the list of possible superconducting phases arising 
just below Tc. 

First of all, the superconducting phase transition corresponds to a particular irreducible rep- 
resentation of the point group G. More exactly, the order parameter AQ,/3(k) should be expanded 
over the basis functions $^^(k) of the given irreducible representation in the vicinity of the critical 
temperature: 




(1) 



The superconducting order parameter is the condensate wave function 




(2) 



(3) 
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Another important property of Ao^(k) follows from the anticommutation relation on fermion op- 
erators in d^): AQ,^(k) = — A^a(— k). Therefore, for crystals with space inversion symmetry, the 
superconducting pairing occurs either in the singlet, space-even states with A(k) = icryipi)^) (where 
ipik) = ■?/;(— k)) or in the triplet, space-odd states with A(k) = icry{a--d(k)) (where d(k) = — d(— k)). 
As the gap in the quasiparticles energy spectrum is propotional to the superconducting order pa- 
rameter, the zeroes of the latter play important role for low temperature thermodynamics of a 
superconductor [^^-ssj ^ ^jj-^g symmetry analysis [®' ''^] shows that in the triplet case the types of zeros 
in the energy gap always correspond to points on the Fermi surface, whereas whole lines of zeros 
are possible for the singlet pairing. 

When this phenomenological approach is applied to heavy fermion superconductors, the following 
one- (ID) and multidimensional (2D or 3D) irreducible representations for different point groups G 
are usually considered (we use notations of p^]): 

ID: Ai, A2, Bi, B2 and 2D: Eig, Eiu, for the hexagonal group Dqu of UPts, UPd2Al3, UNi2Al3. 

ID: Ai, A2, and 2D: E for the tetragonal group D^^ of CeCu2Si2, URu2Si2. 

ID: Ai, A2, 2D: E and 3D: Fi, F2 for the cubic group Oh of UBeig. 

The explicit form of basis functions of these irreducible representations is given in [^^-ssj _ 

Near free energy can be obtained as an expansion in powers of the order parameter which is 
invariant under the action of group Q. 

Consider first the case of ID irreducible representation of G. The symmetry of the order parame- 
ter is determined by the basis function "i>(k). When this representation is identical, ^(k) has the full 
symmetry of the group G x R and corresponds to the conventional Cooper pairing. For nonidentical 
representations corresponding to the unconventional pairing some symmetry elements of G change 
phase of $(k) by tt, and, consequently, the symmetry group of superconducting state is expressed 
in the form G{G') x R, where G' is an invariant subgroup of G From this follows that the 
modulus |$(k)| still has the full crystal symmetry and the superconducting state is nonmagnetic. 
In both cases the order parameters are one-component and their space variations are described by 
one complex amplitude \&(r). For many purposes it is sufficient to write the free energy density up 
to the forth order in powers of ^l/(r). As a result one obtains the usual Ginzburg-Landau functional: 

F = -a{Tm' + P\^\' + K\Dk^\' + ^-^ , (4) 

OTT 47r 

a{T) = a' (^1 - , Dk = dk~i'^Ak, k = x,y,z, 

where a, (3, K are positive phenomenological constants. Magnetic induction B is a spatial average 
of the microscopic magnetic field: B = (h(r)) (h = rot A). Note that for uniaxial (biaxial) crystals, 
the coefficient K should be replaced by a uniaxial (biaxial) second rank tensor. 

Thus the magnetic properties of one-component unconventional superconductors and conven- 
tional superconductors are similar in the framework of GL theory. For this reason we will be 
interested only in multicomponent unconventional superconductors which correspond to the multi- 
dimensional irreducible representations of G. 
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In our review we will consider magnetic properties of unconventional superconductors mostly 
on the example of two-dimensional irreducible representations Ei of group D^h with basis functions 
$3;(k) = k^kz, $j/(k) = kykz in the case of even order parameter Eig, or $a;(k) = zk^, ^y(X) = zky 
for odd representation Eiu- Vector rji in (3) has two components rjr^ and r^y, and the GL functional 
can be written as follows: 

+ KsD*r)*D^r)i + K^Dlr]*D,7], + ____, z, j = x, y , (5) 
A(k,r) = r/,(r)<l,(k)+r7,(r)<lj,(k) . 

Our choice is explained by two reasons. 

(i) The symmetry group D^h of functional (|^) is the same as of the normal state of the heavy fermion 
superconductor UPts. It is the Ei irreducible representation that is involved in many models of 
interpretation of the complicated phase diagram of UPts [24,28,45,74j^ 

(ii) The 2D model (|^) is, in a sense, the simplest example of multicomponent order parameter which 
possesses the nontrivial properties of the mixed state. Many properties of the GL functional (||) may 
be straightforwardly transferred on another more complicated models of superconducting pairing 

in UPts [19'53,83]^ 

The symmetry of the multicomponent order parameter is determined by the coefficients rji in 
expansion (3). Therefore, in order to find the residual symmetry of the superconducting state, 
the corresponding GL functional should be minimized. The complete list of the possible residual 
symmetries for the unconventional superconductivity in heavy fermions is presented in p2.65j ^jj^g 
residual symmetries of the uniform superconducting state in Ei model will be discussed in Sec. 1.3. 



1.2 Magnetic Properties of Usual Superconductors 

One of fundamental properties of superconductors is the expulsion of magnetic field from a sample 
(the Meisner effect). Depending on the material, the Meisner effect can be either complete, as in 
type-I superconductors, or partial, as in type-II superconductors. In the latter case magnetic field 
is dispersed through the material in the form of fiux lines or vortices, each of them containing a 
unit of magnetic fiux 0o = hc/2\e\ = 2.07 ■ 10^^ G-cm^. Note that heavy fermion superconductors 
demonstrate properties of the type-II superconductors. 

Having in mind the discussion of the mixed state in unconventional superconductors, we survey 
the main magnetic properties of usual type-II superconductors with s-pairing, which can be under- 
stood on the basis of the GL functional (^. Varying energy (|^) with respect to S'if* and 5 A, one 
gets the GL equations: 

-a^ + 2/3|^|2*-irAA^ = , (6) 

4-TT 9p Rp'^K 

(V X (V X A)) = — j , 3 = 4k (^*V^ - ^V^*) - 2^ l^l'A . 

c zh he 
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According to the classical results of Abrikosov [^] , magnetic properties of type-II superconductors 
are quite different in different parts of the H-T phase diagram. 

For magnetic fields smaller than the lower critical field Hd, circular surface currents screen the 
external field and magnetic flux does not penetrate inside the sample. The amplitude of supercon- 
ducting order parameter is fixed: |^oP — Q;/2/5, whereas the phase of is arbitrary. The gain of 
the energy in superconducting state defines the thermodynamical critical field He'. Fn — Fg — H^/Stt. 

For H > Hci the flux begins to penetrate inside the material in the form of vortices. A single 
vortex line has an axially symmetric structure: 

*(r) = *(r) exp(-i(/7) , h(r) = h{r) z , (7) 

(r and (p are the polar coordinates) and carries one quanta of magnetic flux [i6-59j ^j^g amphtude 
^(r) is suppressed to zero at the vortex center. It restores to the equihbrium value at a radius 
of the vortex core, which is of the order of the coherence length 

-) ■ (8) 

The magnetic field is maximal at the vortex center and approximately equal to 2i/ci- The super- 
currents (which are beyond the vortex core) screen the magnetic field and h(r) tends to zero outside 
the screening region which is of the order of the London penetration depth 

x4J^T (9) 

The ratio of the penetration depth to the coherence length defines a dimensionless temperature 
independent (near T^) GL parameter k, whose value determines the type of superconductor. For 
typc-II superconductors k > 1/ V^. It is the screening supercurrents jg = c0o/ (Svr^A^r), ^ -C r -C A 
which bring for large k the main logarithmic contribution to the vortex energy per unit length or 
"flux line tension:" 

where the small numerical constant e 0.08 is due to the vortex core energy. The lower critical 
field is defined by the expression: 



^ci = --s. = -^ ln- + 6 . (11) 




When the intervortex distance ri is larger than A, it is sufficient to take into account interaction 
only between nearest neighbors through pair potential U = (j)o/47i-h{ri). In the intermediate region 
^ ^ < A (H > He) vortex lines form a dense lattice. Single vortex interacts with all neighbors 
in the region of the area A^. 

When increasing field reduces the distance between vortices to the order of ^, the superconduc- 
tivity is destroyed. To calculate this the upper critical field Hc2, one should consider the GL equation 
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neglecting the diamagnetic field generated by superconducting currents and the nonlinear term, 
which are small near Hc2- Thus we derive the linearized GL equation: 

= -K {v -i'^^A^ ^, (12) 

which formally coincides with the Schrodinger equation for an electron in a magnetic field. This 
equation is known to have a discrete set of Landau levels. Using the analogy with the Landau level 
quantization, we obtain a set of eigen (critical) fields: 

Tica ^ ^0 1 . . 

" 2\e\K{2n + l) 27i^^ 2n + 1 ' ^ ^ 

The upper critical field corresponds to the highest Hn and, therefore, to the lowest Landau level 
with n = 0: 

- ^ . (14) 

Vortices form a regular lattice of equilateral triangles at the entire field region between Hd and Hc2- 

1.3 New Features of Multi-Component Superconductors 

As we have seen in previous Subsection, the energy functional (^ leads to the particular magnetic 
properties. These properties will be the same (at least near Tc) for all one-component superconduc- 
tors, since they are also described by functional (^) or its anisotropic modifications. Therefore we 
are interested mainly in superconductors described by a multicomponent order parameter. 

The dimensionality of the order parameter is that of the irreducible representation to which 
it belongs. Additional degrees of freedom always correspond to the states of the Cooper pairs 
degenerated with respect to the direction of the internal angular momentum, it does not matter 
whether of spin or orbital origin. For example, in the case of the 2D order parameter, the degeneracy 
exists between the states with "up" or "down" direction of the angular momentum i{'q*xr]). 

In the case of GL energy (|^), the equilibrium order parameter satisfies the following nonlinear 
equations: 

ar,i - 2Pir]iiri*ri^) - 2P2vKVjVj) + K,D,Djr]i + K^D.D.r]^ + KsD,D,r]j + K,D,D,r]i = 0, (15) 

4e r ] 
ji = Im [KirijO^rij + K2ri*Djrij + KsrijDjT]^ + K46i,r]*D^r]j . 

New properties in comparison with ordinary superconductors appear already without magnetic 
field. Depending on the sign of (32, the resulting superconducting phases are: 

/ \ 1/2 

A(k)=f^j (^i(k)+2$2(k)), (|4i| = |$2|,^l±^2), /32>0, (16) 
/ \ 1/2 

^^^^=(2^;) (* = **), P2<0. (17) 
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To be consistent, the energy should be bounded from below for each phase (|T6|) and ( |T7|) what leads 
to the positive definiteness of the forth order terms: 



/3i > , /3i + /32 = A2 > . (18) 

For /?2 > the ground state of the GL functional has two-fold degeneracy with respect to 
the appearance of either $_|_ = $1 + i$2 or = <&i — phases. Both phases are conjugated 
with respect to the time inversion. Physically, these phases correspond to the superconducting 
states with definite values of the projection of internal angular momentum of Cooper pairs like 
in ^He-A. Symmetry group of $±(k) phase includes arbitrary rotations about z combined with 
multiplication by the phase factor exp(=pz9?)L^, and rotations by an angle of vr about x combined 
with the time reversal RU2x- Discrete degeneracy of the energy minimum leads to the two types of 
superconducting domains. 

For /?2 < the order parameter is real up to the complex factor. The symmetry group of the 
superconducting state in this case is D2{C2) x R. The direction of the vector $(k) is fixed by the 
sixth order term ri*^rj^ +7]*^!]^ in the GL functional (?7± = rjx^irjy). We have not included this term 
in d^) since it is small near in comparison with the other terms. In this case the ground state 
possesses continuous degeneracy with respect to the direction of Therefore smooth texture of 
the order parameter A(k, r) = n(r) ■ $(k) should be expected in real samples rather than domains. 

Derivation of the free energy functional (|^) in the framework of the Landau theory requires as a 
necessary condition positiveness of the gradient terms without any magnetic field (i.e., if A(r) = 0), 
which leads to: 

i^i > , > , 7^123 = K^ + K2 + K^>Q . (19) 

However within this range of parameters it is possible the magnetic instability of the uniform 
superconducting state p^], which results in a spatially varying order parameter in zero applied field 
[®^]. In fact exact restrictions connected to the existence of the lower energy boundary for such 
laminar phases do not differ strongly from the above necessary conditions p^]. Therefore in the 
following consideration we impose only conditions (0). It will be useful to introduce dimensionless 
parameters C = {K2 + K^)/2Ki, D = {K2~K^)/2Ki. Then, according to (|l|), we have 1 + 2^ > 0. 

Neglecting by influence of the order parameter A on the pairing potential, one can obtain in the 
weak-coupling approximation the following relations between phenomenological parameters for Eig 
representation [^^]: /?2 = 0.5/?i, Ki = K2 = K^. Note that the approximate particle-hole symmetry 
of Bogoljubov equations near the Fermi surface results in a relatively small value of D ~ {T^/ EpY 

|-22,14,55J 

Returning to the magnetic properties of unconventional superconductors, we could expect the 
same qualitative features in their behavior, namely, penetration of separate vortices inside the 
material at fields higher than Hd-, formation of the vortex lattice, which becomes denser with 
increasing external field, and destruction of superconductivity dX H = Hc2 when the distance 
between vortices becomes of the order of ^. 
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However, the GL equations (|1^) for the 2D order parameter have a more comphcated structure 
than those for the s-wave superconductor (^. This is why new properties of the mixed state of 
unconventional superconductor are expected, and the well established results for usual supercon- 
ductors should be reexamined in the multicomponent case. 

The structure of the mixed state of a multicomponent superconductor is considered in the 
following Sections of this review. Here we give a short outline of the problems to be discussed. 

1. The multicomponent structure of the order parameter and two different types of ground state 
lead to the possibility of the existence of different types of vortices carrying an integer flux, which 
are discussed in Sec. 2. Unlike the case of usual superconductors, the vortices can have a nonsingular 
structure of the core and be nonaxisymmetric. Since in reality only the most energetically favorable 
vortices exist, one should calculate the energy for each type of vortices and choose the vortex 
with the lowest energy. Note that the favorability of single-quantum vortices over multiple-quanta 
vortices is also questionable. 

2. In order to calculate the upper critical field Hc2(T) for unconventional superconductors, the 
multicomponent linearized GL equations have to be solved. The puzzling feature of solutions of 
such equations described in Sec. 3 is that all components of the order parameter participate in the 
formation of the Abrikosov vortex lattice, and, hence, the internal structure of the mixed state near 
Hc2 is rather complicated. The main question to be considered in this context in Sec. 4 is: whether 
the lattice has the regular triangular form as for s-wave superconductors or it is distorted. Another 
surprising fact is the existence of sevral superconducting phases from different Landau levels with 
close values of critical fields. It will be shown in Sec. 5 that the admixture of such phases inevitably 
leads to structural phase transitions in the Abrikosov lattice below Hc2- 

In addition to the review of already published works we have included new original results in 
Sees. 3.5 and 4.3. 

In fact, depending on the parameters in multicomponent GL equations, both nontrivial vortices 
near Hd and structural phase transitions in the Abrikosov lattice near Hc2 may occur. The crucial 
open question of the mixed state theory of the multicomponent superconductors is relation of these 
two limiting cases: 

— How are nonaxisymmetric vortices placed in the lattice near -ffci? 

— What is the behavior of such vortex lattice with the increase of magnetic field? 

— What kinds of phase transitions can occur in the intermediate region between Hd and Hc2 in 
unconventional superconductors? 

2 Vortices in Unconventional Superconductors 

The vortices in a superfiuid liquid with the multicomponent order parameter were intensively studied 
during the last two decades in connection with rotating superfiuid ^He (for review see It was 
shown that the response of such a liquid to rotation is formation of topologically stable vortices of 
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the order parameter. At large distances from the vortex core the slow space variation of the order 
parameter A(k, r) can be described as a long-range texture of the uniform order parameter A(k) 
which is fixed by the forth order terms in the GL functional. (For ^He it is the space variation 
of A or B phase.) However, in the vicinity of the vortex core the bulk phase is destroyed and all 
components of the order parameter have the same order of magnitude. We define the characteristic 
size of this region as ^. 

The analog of the superfluid liquid rotation is the response of the superconductor to the external 
magnetic field. Before translating the results of superfluid ^He theory on the case of the multicom- 
ponent superconductors, we should stress that there is a feature which makes the difference between 
superconducting and superfluid vortices: the first ones have the finite size A due to the screening 
effect of the currents, whereas the size of superfluid vortices is infinite. 

Therefore, three length scales enter the problem of vortices in multicomponent superconductors: 
the vortex size (the London penetration depth) A, the radius of the vortex "hard" core (the coherence 
length) ^, and the radius of the region around the vortex core ^ where all k-components of the order 
parameter are mixed. Note, however, that such a stratification of the vortex on the characteristic 
scales has physical sense only under the condition that A ^ ^ where the values A, C,, C, are 
defined by the parameters of multicomponent GL functional. Having in mind this case, we will 
discuss magnetic vortices in the two-component model (^, which have been studied in most details. 
We will assume also that the magnetic field is parallel to the principal crystal axis. 

The properties of vortices and the scales of A, ^, ^ strongly depend on the bulk superconducting 
state. As we have already mentioned in Sec. 1.3 there are two stable minimums of the GL functional 
(^ describing homogeneous state for different signs of /?2- There is either a complex order parameter 
(broken time-reversal phase) A(k) ~ $i + i^2, {P2 > 0) or a real one (vector phase) A(k) ~ 
{P2 < 0). We will consider both cases separately in Sees. 2.1 and 2.2. 

2.1 Vortices in the Broken Time- Reversal Phase 

To begin with the discussion of vortices in the broken time-reversal phase, note that the ground state 
is doubly degenerate with respect to the direction of the angular momentum of Cooper pairs. Due 
to this discrete degeneracy, the domain structure (with phases $+(k), $_(k)) in the ground state of 
the two-component superconductor is expected . Thus, the vortices originating from both phases 
can exist inside the superconductor. Note also that the domain structure itself contributes to the 
magnetic properties of the unconventional superconductors because of two important effects: (i) the 
internal ( "up" and "down" ) magnetization of Cooper pairs owning to their orbital momentum (this 
bulk magnetization is assumed to be small due the microscopic electron-hole symmetry [''^] ; (ii) the 
magnetization owning to persistent surface currents which circulate around domains and originate 
due to the inhomogeneity of the order parameter (on the scale ^) near the domain boundary (these 
currents are due to boundary effects of the bulk magnetization [''^]). The low field magnetic response 
of unconventional superconductors due to the domain structure was considered in detail in p3,64,65j 
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Note that an external magnetic field lifts the $±(k) degeneracy and makes one sort of domains 
more favorable. Therefore, the ground (vortex) state of unconventional superconductor in an exter- 
nal field is rather complicated and history (field or zero- field cooling experiments) dependent. In 
this review we will neglect for simplicity the effects of domain- and internal-magnetization and will 
assume that the vortices exist in the bulk of domains (the latter means that the domain sizes are 
larger than the penetration length). We will discuss the bulk vortices in both phases $+(k) and 
^_(k) simultaneously using "±" subscript. 

2.1.1 Global vortex structure 

The long range structure of n-quantum vortex is given by: 

A±(k, r) = M/±(r)($i(k) ± z$2(k)) , (20) 

where phase of \l'±(r) varies on 27m about the vortex center. Substituting (^) into (0), we obtain 
the effective GL equation for the space variation of \l/±(r): 

- a ^ ^ h{r)DKi) ^± + 4/3i|^±|2^± - Ki{l + C)DiDi^!^ = . (21) 
nc J 

The term proportional to D is responsible for the above mentioned interaction of the angular 
momentum of Cooper pairs with local magnetic field. If we neglect by this term (see Sec. 1.3), Eq. 
(^Tj) becomes the same as the usual GL equation (|) for s-wave superconductor. Then by analogy 
we find the size of the vortex (penetration depth) in the time-reversal breaking phase: 

ll6jreWs:i(l + C)i ' ^ ' 



We can also define the coherence length which corresponds to the size of the vortex core: 



1/2 



However, this length has a limited physical sense for small /?2 because in the region r < other 
components of the order parameter are admixed, and Eq. (^TJ) does not hold. In this region the 
gradient energy due to the texture of $± 



K,{l + C) 

^grad ~ ^bulk [^^) 

St 



becomes comparable to the contribution to the uniform part of the GL energy due to other admixed 
components of the order parameter: 

SFbuik ~ ^ Fiyuik ■ (25) 
Pi 
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Equating SF^y^ik to Fgrad, we obtain the characteristic scale ^t'- 

(26) 

The same estimation is vahd also for the width of the domain wall between #_|_ and phases 

By the analogy with usual s-wave superconductor, we can conclude that the one-quantum vortex 
has the smallest energy and Hd is defined by formulas (|^) and ( [Tl| ) with a formal substitution ^ C,t, 
e —>■ if, where it is the energy of the vortex core in the region r < . In the region r ^ the vortex 
(pOD is described by the same function (0): \E'±(r) = \l'(r)e~*'^ (with substitution ^t, ^t) like a usual 
Abrikosov vortex. 

These statements are valid provided ^ ^t, which is the case for not very small positive ^2- In 
particular, we have assumed that the main vortex energy is concentrated in the screening currents 
region At > r > ^t, what means that ln(At/^t) ^ if. 



2.1.2 Vortex core structure 

The investigation of the vortex core structure at > was done in [™''^i>6j_ Our description of the 
vortex structure is based on the results of these papers. We will use the analytical arguments of 
to illustrate the appearance of vortices with nontrivial core structure and then present the results 
of the numerical calculations of [™] which were done in a broad range of parameters. 
In the region r < C,t all components of the order parameter are mixed: 

A±(k,r) = ($i(k) ±2$2(k))*±(r) + ($i(k) ^ z$2(k))X^(r) . (27) 



Functions \E'-i-(r), X^{r) can be determined from solution of the GL equations (|T^). Since at r 3> 



It 



vortex ( P7| ) should gradually turn into (^^, we have at the large distances \E'±(r) — > \E'(r)e and 
X^(r) ^ 0. 

Axisymmetric vortices 

In general, there exist different solutions of Eqs. ([15|) obeying the boundary conditions \E'±|r._^oo = 
\J a/A(3i e~*'^, X^\r^oo = 0. One of the extremums of GL functional always corresponds to the "most 
symmetric vortex, " which possesses all symmetries possible for line defects in condensed matter [^°]. 
According to [®°] the maximal symmetry group of a line defect includes continuous rotations about 
its axis with multiplications by a phase factor, the reflection ah in the perpendicular plane, and the 
combined symmetry RU2- The last two discrete symmetries have no effect on vortex classification 
in the model of 2D superconducting order parameter. 

The generator of continuous subgroup and corresponding condition of invariance have form: 

gA^(k,r) = 0, Q = L,-qi. (28) 
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Here / is the generator of U{1) group, and q is an integer which has the meaning of the total 
angular momentum eigenvalue of the vortex. The symmetry group of the axisymmetric vortex is 
{e~^'^'^L^, RU2, d'h}- The solutions of Eq. ( P^ for the case of one-quantum vortices in the phases 
^± are 

vl>- (r ) = v]>+ (r ) e-^ , (r) = X_ (r ) e*'^ , g = , 

^^{r) = ^_{r)e''^ , X^(r) = X+(r) e'^*^ , q = -2 . 

Functions \E'±(r) and Xzp(r) vanish at r = making a "hard" core of the axisymmetric one-quantum 



vortex. At large distances I^E'il ~ a/A/3i, and amplitude X:^ is small. Linearizing Eq. (|T5D in X: 



in this region, one comes to the following equation: 

( - « + 4(A + 2/?2)|^±|' - K,{1 + C){Dl + DD) = K.CiD, ± zD,)2^± . (30) 



Substituting = ya/4/3i e in (pO]), one obtains: 

X^r) ^ - 




From ( pTD one can easily estimate the region where X" becomes comparable to For C > 1 its 
sizes coincide with In order to find behavior of \l'±(r) and X=p(r) at r ^ 0, one should solve 
both GL equations (|T3p simultaneously. The admixture of the other components changes the usual 
assimptotic '^±{r) ~ r of superconducting amplitude at the vortex core. 

Noricixisymmetric vortices 

Although the axially symmetric vortex solution is always extremum, it is not necessarily the 
absolute minimum of the GL functional. Obviously, it is the minimum as .^t — {j32 ^ Pi), 
when \E'^(r) is defined by Eq. (|30|) almost throughout the entire space. Functions \l'^(r) are so- 
lutions of (|30D in all space also for C = when the additional amplitude in axisymmetric vortex 
X^(r) = 0. Due to this feature the case of C = is a convinient starting point for consideration of 
nonaxisymmetric vortices. 

In this case Eq. ( |30D becomes uniform and takes form of the eigenvalue problem for the ordinary 
Schrodinger equation on Xzp in a potential well with f/(r) = —a + 4(/3i + 2/52)\E'^(r). A two- 
dimensional well always possesses a bound state, which is placed on some finite energy below the 
value of potential at infinity U{r)\r-,oo = a/52/2/3i > 0. Therefore at C = the lowest energy level 



of (|30D should pass through as decreases. This means that the vortex A^(k, r) = \E'j.(r)^>-|-(k) 
becomes unstable towards perturbation of the form X°(r)$^(k). Here X^{r) is a real axisymmetric 
function exponentially decreasing at r ^ and with maximum at r = 0, which corresponds to the 
lowest level in the well. 

The critical value of f3 = f32/(3i corresponding to such an instability in the vortex core is f3c = 0.24 
according to [^^] . Using trial functions Barash and Mel'nikov [^] have found a slightly different value: 
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0.37. However, subsequent numerical investigations of one-component odd-parity model [^^], which 
reduces in a magnetic field to model @ with C = 0, have suggested the same value for Pc as in 
The order parameter amplitude at the center of such nonsingular vortex grows as (/?c ~ PY^"^- 

When C is exactly equal to zero, sum |\E'-|-(r)p + |X°(r)p depends only on r. This leads to axial 
symmetry of the whole vortex. Such a symmetry of the nonsingular vortex, however, is "accidental" 
and disappears for any small value C. 

For small C solution of Eq. (|3D|) is the superposition X^(r) = X"(r) + X°(r). Consequently, 
X^(r) is not already invariant under arbitrary rotations and axisymmetry is spontaneously broken. 
Moreover, using (|29D we come to the following conclusions about symmetry of the nonsingular 
vortices: 

(i) Vortices originating in state are one-quantum, axisymmetric and singular for P2/ Pi > Pci^)- 
For P2/P1 < Pci^) they have nonaxisymmetric triangular stiuctuie (Fig. 2a). The residual sym- 
metry group is: {e'^'^'^^^L2n/3, <5"/i, Rd'y}, where (j„ means the refiections with respect to the triangular 
medians. 

(ii) One-quantum vortices in state have axisymmetric singular structure for P2/P1 > Pc{C). For 
1^2 / Pi < Pti^) vortex cores are of nonaxisymmetric "crescent" shape (see Fig. 2b). The residual 
vortex symmetry is of the vector type and contains only two elements: {ah,Rcry}, where is the 
reflection in the symmetry plane of "crescent." 

The boundary of the axisymmetric vortex instability Pf{C) was found by Tokuyasu, Hess and 
Sauls [^°] numerically (see Fig. 3). The results of at small C qualitatively coincide with [™]. 

For topological reasons, only the vortices with an integer number of flux quanta can exist in the 
bulk of time-reversal breaking phase. However, vortices carrying fractional number of magnetic flux 
quanta are known on a domain boundary between phases and p^]. Moreover, two fractional 
vortices could have a lower energy than a single-quantum vortex. Starting on this picture, Izumov 
and Laptev p^] suggested a scenario in which a bulk vortex decays forming a small domain of the 
time-conjugated phase and two fractional vortices which exist on the closed boundary of the new 
domain. The loss in the gradient energy due to the domain boundary is compensated by the gain 
of energy due to the vortex decay. Note that the vortex core instability in the phase ^+ with the 
formation of the nonaxisymmetric crescent structure and phase inside the vortex core is an 
extreme case of the Izumov and Laptev scenario when the domain size is sufficiently small. 

Developing the analogy with vortices in s-wave superconductor established by Eq. (21) we can 
estimate vortex core energy It as {\n.^t/it+^)- Then for those P2, which satisfy condition ^t/it > ^t/it 
or P2 < P\l \fK^ the main part of the vortex line tension is the energy of the vortex core q rather 
than the hydrodynamics part proportional to ln(Af/^t)- In this case the energetic favorability of 
the one-quantum vortices in comparison with multi-quantum ones is broken. For example, the 
two- quantum axisymmetric vortex in the phase has according to (pSf) nonvanishing amplitude 
X'^{r) at the center and therefore smaller core energy it than one-quantum vortex. Its existence in 
the particular parameter region was found numerically in [''^]. 
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Spontaneous breaking of vortex axisymmetry was known before in superfluid ^He-B. In particular 
parameter region the core of so called v vortex is subjected to the '^director type" instability f'^]. As it 
have been shown above the axial symmetry breaking in the vortex cores is quite general phenomenon 
in the case of the multicomponent order parameter. Ordering of such nonaxisymmetric vortices in 
the lattice at H > Hd is still an open question. 

2.2 Vortices in the Vector Phase 
2.2.1 Global vortex structure 

For /?2 < 0, the equilibrium order parameter is A(k) = (n ■ $(k))\E', where n is a unit vector. Far 
from the vortex core the long range texture of the order parameter preserves this form: 

A(k,r) = (n(r) ■$(k))^(r) . (32) 

Substituting ( p^ ) into (||) and keeping the gradient terms responsible for Goldstone modes, we have: 

F = -a\^\^ + pu\'^\'' + KiD*^*Di^ + K23{niD*^*){njDj^) + 

+ |^P(iri(rotn)2 + iri23(divn)2) , (33) 

There are two topological reasons for the formation of stable line defects in the vector state (0): 
(i) the phase of \I'(r) changes on 27m in making a complete circuit; (ii) the vector n(r) sweeps the 
unit circle m times when moving around the vortex. Therefore, possible defects are classified by two 
quantum numbers {n,m), both which are integer or half-integer (for general aspects of topological 
classification see [^^]). 

Note, however, that the phase variation about the vortex leads to supercurrents which screen 
the magnetic field outside the region of the size A. For this reason the vortex energy has the 
upper cutoff parameter A and is therefore finite. The variation of the vector n(r) does not produce 
screening currents. Therefore, for m 7^ 0, the contribution to the defect energy ~ / (i^i(rot n)^ + 
Ki23{dw n)'^)dV is logarithmically divergent beyond the core. It makes the defects with m 7^ less 
favorable then the pure phase vortices with m = but n ^ 0. 

Since we are interested in the vortices with minimal energy, we will consider only those with 
uniform distribution of n(r) at large distances. Assuming that n(r) is uniform, we rewrite the 
effective GL equation for A(k, r) = \E'(r)(n ■ $(k)) as 

- + 2/?i2|^|^^ - KiDiDi^ - 2KiC{niDif^ = . (34) 

This equation has a form of the one-component GL equation for a crystal with anisotropy axis n 
and magnetic field H directed perpendicular to n. Using scaling transformation of p^]: 

ri, = r(, v^l + 2C , r± = , (35) 
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we reduce (p^) to the isotropic GL equation 

-a^ + 2/5i2|^|^^ - KWl + 2CD[D[m = . (36) 
From (^) we conclude the following: 

(i) the one-quantum vortex {n = 1) is the most favorable one; 

(ii) the sizes of the vortex corresponding to the anisotropic London penetration length are: 

A. 



where 



All = A„ v^l + 2C , A^ = " , (37) 



(iii) the coherence length corresponding to Eq. ( p6D is defined as 

1/2 



a 



(iv) at distances larger than ^„ the vortex has the form: 

A(k, r) = ^(r)(n ■ $(k)) = ^(r')e-^'^' (n ■ $(k)) , (40) 

where function \E'(r') has the same dependence as the usual Abrikosov vortex (|^. As in the case 
of the vortices in the broken time-reversal phase, the characteristic scale of is calculated by 
comparing the contribution to the free energy due to admixing of the ^± component and the 
characteristic gradient energy on the scale ^t,. Since the density of the gradient energy is anisotropic, 
the region where all components of the order parameter are mixed is defined by the anisotropic scales 

U-L</TT2C, a = (41) 

where 

The lower critical field for the vortices in the vector phase is calculated from ( pT]) by substitution 
A — s> At,, ^ — * ^t,, e — > e^. Note again that the above description of the long range vortex structure 
is valid only if A^, ^ 

2.2.2 Vortex core structure 

The calculation of the vortex core structure for /?2 < has not been done yet. Therefore we will 
outline only expected qualitative features of them. 

Similarly to the case of /32 > 0, we start from the description of the "most symmetric vortex." 
In order to find its "full symmetry," we note that the rotational symmetry of the vortex is broken 
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due to the vector n: {L^^, RU2X, P}- The general structure of the "most symmetric vortex" inside 
the region is given by 

A(k, r) = (n • *(k))V'(r) + z(m • *(k))x(r) , (43) 

where m denotes the unit vector which is perpendicular to n and H. Functions ip{r), x(r) have D2 
symmetry, and their phases vary by 27r about r = 0. Outside the core region ^(r) \E'(r')e^*'^', 
x(r) 0. Both amplitudes ■0(r) and x{'^) vanish at the vortex center. 

If C|h ^-L order of {y, the "most symmetric vortex" is the most favorable. When 

^„ becomes sufficiently large in comparison with one can expect the appearance of vortices with 
a lower symmetry as is the case with the vortices in the broken time-reversal phase. However, this 
question requires further investigations. 

2.3 Discussion 

Summarizing this Section, let us discuss what conclusions we can draw from the outlined picture of 
vortices in unconventional superconductors and what problems they pose for further studies. 

The appearance of the unconventional vortices was mostly studied in the framework of the 
two-component model when H || c. It was shown that in some range of parameters the vortices 
originating from the broken time-reversal phase {j32 > 0) have the nonaxially symmetric structure 
of the core, unlike the vortices in usual superconductors which are always axisymmetric. It was also 
shown that the domain structure of the ground state of a superconductor plays an important role 
in its magnetic properties. 

The instability of the vortex core with axisymmetry breaking, the domain formation and the 
vortex interaction with domain boundaries are certainly the features connected to each other. How- 
ever, the most favorable ground state of such a system in a finite magnetic field near Hd and, in 
particular, the exact symmetry of the lattice of nonaxial vortices are still unknown. 

The vortices in a vector phase {f32 < 0) were studied much less. It is known that the vortices have 
an anisotropic shape in the direction of ^(k). However, the structure of the vortex core (depending 
on the parameter (32) needs further investigations. 

In the foregoing the vortices have been treated in the limit when a vortex can be clearly divided 
into two principal regions: the region of the screening currents, in which the main vortex energy is 
concentrated, and the vortex core region, in which the axisymmetry instabilities occur. This limit 
corresponds to the restriction A ^ ,^ which is the case whenever (32 is not too close to zero. Note, 
however, that this case is in a sense the most trivial one because all the results were obtained by the 
straightforward generalization of two famous theories: Abrikosov theory of the vortex structure (for 
the region r > ^) and the theory of the vortex core instability in the superfiuid ^He [^°] (adopted 
for the two-component order parameter $1 ± or 

A very interesting problem comes out from the above discussion: what is the structure of the 
vortices at small P2 when the condition A > ^ is not satisfied? In this case not only the currents 
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which originate from the gradient of phase of the order parameter but also the currents due to 
the variation of all other admixed components of the order parameter contribute to the screening 
of magnetic field. Another related question is how the vortices in the broken time-reversal phase 
transform to those in the real (vector) phase when the parameter changes from > to /32 < 0. 
Note that these problems require a careful consideration because even the ground state of a two- 
component superconductor is nonuniform at small (52 {H — 0) and has a space-modulated laminar 
structure [^^'^^]. Generating of vortices in this state (and therefore the problem of Hd) is also not 
clear enough p9'56j 

As we have seen, magnetic vortices in a two-component superconductor have a lot of unusual 
features. However, the two-component model is the simplest one which describes the unconventional 
superconductivity. Therefore, one can expect even more peculiar new features in other multicompo- 
nent superconductors. In this connection, it is worthwhile to mention the mixed state in the triplet 
superconductor with a weak spin-orbital coupling which was studied by Burlachkov and Kopnin [^^]. 
They showed that in magnetic field spin-vectors of such a superconductor form a distinct texture 
which allows the applied field to penetrate inside a superconductor without forming singularities. It 
is remarkable that the lower critical field for such a texture can be much smaller than corresponding 
field for the usual singular vortices. 

3 The Upper Critical Field in Unconventional 
Superconductors 

The calculation of the upper critical field for the superconductors with multicomponcnt order pa- 
rameter is analogous to the case of s-wave superconductors (see Sec. 1.2). Starting from the mul- 
ticomponcnt GL energy functional, we derive the corresponding linearized GL equations under 
assumption that magnetic field inside a superconductor is uniform and equals to the external ap- 
plied field. Therefore, the calculation of the upper critical field reduces to the search for eigenvalues 
of the system of the linearized GL equations and then to the choice of the critical (eigen) field with 
the maximal value. 

In Sec. 3.1 we present characteristic examples for which linearized GL equations for multi- 
component order parameter admit analytical solutions. We will see that, despite the complicated 
structure, their eigensolutions can be classified by definite quantum numbers. The nature of this 
classification is a consequence of the symmetry of the normal metal in magnetic field (Sec. 3.2). 
Hence, properties of the superconducting phases can be discussed without solving the linearized GL 
equations and general conditions of a kink on the Hc2{T) curve (the problem related to the phase 
diagram of UPts) can be formulated (Sec. 3.3). In Sec. 3.4 we will enumerate the specific features of 
the Hc2 anisotropy which make it possible to distinguish the unconventional superconductivity from 
the usual one. Finally, in Sec. 3.5 we consider the possibility of the appearance of superconducting 
phases modulated along the field direction. 
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3.1 Upper Critical Field in the Two- Component Model 

In this Subsection we consider two cases with H || c and H ± c when hnearized GL equations 
for the two-component order parameter can be solved analytically [ii'67,74,79j_ rjj^g question of 
anisotropy for arbitrary directed field will be considered in Sec. 3.4. 

To obtain the linearized GL equation, we neglect in (JT^) the forth order terms. For H || z these 
equations are 

I'fJ I C{D,-iD,f (i + C){Dl + Dl) + f^HD)W)' ' 

where ri± = rj^ i: irjy. In ( ^if ) we set DzT] = 0, since the variation of the order parameter along H 
only decreases Hc2- Defining the creation and annihilation operators 

1/2 

{D^ - iDy) , \a, a+j = 1 , (45) 



he 



one obtains the eigenf unctions of (^41) as combination of the r-dependent Landau levels wave func- 
tions /n(r) (for which a^a fn = {n + \)fn) and k-dependent basis functions $j(k): 

A„(k,r) = /„(r)($,(k)-^<l>,(k)) + /„_2(r)(l'.(k)+i$,(k)) , 



A„ = (2n - 1)(1 + C) ± ^(2 + 2C - Df + AC^n{n - 1) . (46) 

Their respective eigenfields if„ are related to the eigenvalues A„ of the linear problem by the 
expression 

Hn = , (47) 

and Hc2 = max{if„}. The solution with the lowest eigenvalue is either 

A„(k,r) = /o(r)(|..(k)-^<|.,(k)) , 

A, = Ao = 1 + C - D , for D > , (48) 



or 



SKiKr) = /o(r)(<l,(k)+z<l,(k))+^/2(r)(<l>,,(k)-z<l>,(k)) , 



XsK = A2 = 3(1 + C)- ^8C2 + (2 + 2C - D)2 , 
u;= ^'~^~^~^ , foTD<-^. (49) 

2V2C 1 + C ^ ^ 



The new feature of the anisotropic pairing which can be seeing from ( ^61 ) is the nonfactorized 
dependence of A on k and r. Scharnberg and Klemm (SK) were the first who discovered this fact 
solving Gor'kov equations for p-wave pairing in the weak-coupling limit [^^]. These solutions were 
obtained within the phenomenological Ginzburg-Landau approach in [6''>79j Note that the solutions 
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( ^q) contain those products of k- and r-dependent functions, for which the sum of the projection of 
the Cooper pair angular momentum m and Landau level number n is the same. The corresponding 
quantum number N = n + m was introduced in [^^] . We shall call N as "generalized Landau level 
number." The symmetry reasons for the appearance of such quantum number will be given in the 
next Subsection. 

For different values of phenomenological constants C, D, the smallest eigenvalue is either Xa or 
XsK- If the approximate particle-hole symmetry exists near the Fermi surface (e.g. in the weak- 
coupling regime) then D ^ 0. Therefore it is the SK- phase (^) which should appear at ifc2 in the 
"real" two-component superconductor (see Fig. 4). 

Due to the cylindrical symmetry of the GL functional (^), Hc2 is isotropic in the basal plane. 
Taking H || x, we obtain from (|15]) the following system of equations {K = K4/K1): 





2C)Dl + Dl + KDl 2CD,Dy 

2CD^Dy Dl + {1 + 2C)Dl + KD] 




(50) 



This system can be solved analytically only under the condition D^ri = (see Sec. 3.5) when 
two equations in (|50D become decoupled and coincide with the one-component GL equations for 
anisotropic s-wave superconductors p^]. Using the scaling transformation, one obtains in this case 
the eigenfunctions of (|50D with the lowest eigenvalues [^^]: 



Ai(k,r) = fois,y,si'z)^^ik) , Xi = VK , (51) 
.1 = ^, 



A2(k, r) = fo{s,y, s^h) %{k) , A2 = jKil + 2C) , (52) 



S2 = ^K/il + 2C) . 

Therefore, for C > 0, the phase (^) corresponds to the upper critical field, whereas for C < the 
phase (B^) appears in the vicinity of Hc2- 



Scaled zero Landau level functions (|5T1) and (|52D are superpositions of all even Landau levels. 
Therefore the above classification by quantum number reduces to N over mod 2. In addition, 
solutions (1^), (0) and (^) are characterized by the parity with respect to the reflection in the 
plane perpendicular to the filed direction. 



3.2 Symmetry Classification of Superconducting Phases near Hc2 

The appearance of the Abrikosov vortex lattice in type-II superconductors at Hc2{T) is a second 
order phase transition. Hence, it can be considered in the framework of the Landau phase transition 
theory p^] irrespectively of particular model assumptions, such as the GL approximation. The 
corresponding symmetry approach was developed recently in the works [3^,40,821 j^g^g g^jg^ ^^])- 

The appearance of the vortex lattice results in spontaneous breaking of the symmetry group 
of the normal state in a magnetic field. Let us first consider the Cooper pairing in an isotropic 
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metal, when G = SO^ x P. In the presence of the magnetic field H directed along z axis the total 
symmetry group Q is reduced to 



g = Tx D^hiCooh) X f/(l) . (53) 

The group Dooh{Cooh) contains the subgroup Cooh of arbitrary rotations about H with the refiection 
in the perpendicular plane ah and the combined symmetry elements RU2 and Ray. Here U2 are 
rotations through angle vr about axes perpendicular to H, and a^, are reflections in planes containing 
H. 

The phase transition of a continuous type implies that the order parameter obeys the linear 
equation on the transition line: 

C{H,T} A{k,r) = 0. (54) 

This equation can be considered as a symbolic form of either differential GL equations, e.g. ( PD and 
(pOD , or integral Gor'kov equations, see p8'62j_ ^]^q linear operator C is invariant under the action 
of the symmetry group Q. In general, the analytic solution of (|5^) is a complicated problem. But it 
is known that different solutions of linear equations belong to different irreducible representations 
of the corresponding symmetry group. This classiflcation is model independent and valid on the 
whole line Hc2{T). 

Note that the group ( ^3l) has the same structure as group except for the breaking of time 
reversal symmetry. For instance, subgroup T is an abelian group of three-dimensional translations. 
However, in magnetic field symmetry operators acting on the order parameter include gauge trans- 
formations (along with the changes of space variables), which return the chosen vector potential 
A(r) to its original form: 



= exp |-i ^ ^ SkpAp{sqjr'q) - Ak{r'j) dr'^ ■ exp{i(p(i^ + Ik)) , (55) 



where 

4 = -i{xdy - yd^) , 4 = M 7^ X k )^ , Skp = I ^^^"^ smv? \ 

\ ak J \~ simp cos if J 

We wrote the gauge transformation for the Cooper pair as that for a single particle with a total 

charge of 2e only for the sake of brevity. The exact transformation operators acting on the order 

parameter (Q) include exponents of the sum of two integrals, both in the above form, with upper 

limits ri and r2 and multiplied by charge e. Our simplification does not change commutation rules 

for Ta. Therefore, symmetry classification given below is valid not only in the GL regime, in which 

r varies over distances ^ ^ but also at T — when the space variations are of the order the 

Cooper pair size. 

For operators ( [55| ) the following relations hold: 
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TJ'h = ThT^ exp <{ z — H [a X b] 



2e 



(57) 

^•fii+!fi2 ■ (58) 
Although the two-dimensional translations in x-y plane commute with each other, this is not true 



for the corresponding magnetic operators (|57|) . As follows from (p6D, the operators Ta form not an 
ordinary vector representation of the abelian group T, but a ray one (e.g., see [^'^]). Different choices 
of the gauge for the same distribution of the magnetic field correspond to the equivalent ray rep- 
resentations. Near Hc2 the superconducting order parameter is characterized by its transformation 
properties under the action of magnetic operators or, equivalently, by the corresponding irreducible 
ray representation. 

In order to classify the different irreducible ray representations, we introduce the generators of 
infinitesimal ray transformations: 



2e 

tx = -idx - — {A^ + Hy) 
nc 

2e 

ty = -idy- — {Ay- Hx) 

Lz = Iz ~\~ Ik ~ 1~ I xA 

nc 



(59) 



yA, - y (a;' + y^] 



which satisfy the following commutation relations: 

2e 



he 



Hi 



^xi L/z 



—it. 



y ' 



(60) 



In the case of axially symmetric gauge A = (—^Hy, ^Hx,0) magnetic translations are reduced to 
the exponent functions from usual magnetic generators [^''^^]: 



nc 



-idy + ^ X 
nc 



(61) 



while is the ordinary rotation. 

By virtue of (|60D , we introduce the Casimir operator J', which commutes with operators ( ^51) 
and operator C. Its eigenvalues enumerate irreducible representations for a given ray representation: 



J 



4:\e\H 



{ti+t'y)+Lz 



One can find the proof of the commutation relations of and derivation of 
lowering and raising operators (|^), we can write JT" in the another form: 



(62) 



in p ]. Using 



jf = a+a + 4 . (63) 
It is useful to construct the general type of eigenfunctions of operator (and therefore £): 

(64) 



A^(k,r) 



Y: A^{H,T)Ur)^i{k) 

N=n+m 
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where ^^C^) wave function corresponding to the definite projection m of the internal angular 

momentum on z. From the symmetry point of view, each product /ri(r)\i/^(k) with n + m = N 
can be considered as one from the equivalent choices for the basis functions of the A^th irreducible 
ray representation. Then constants An{H,T) depend on the construction of the operator C in the 
particular model. In addition to N, we define another quantum number o" = ±, which is the parity 
under refiection ah in the plane perpendicular to H. 

It can be easily seen that phases (|46|) present a specific case of ([641). The quantum numbers 
(A^, cr), which were introduced in the previous Subsection, are conditioned by the operators J', ah- 
Note that the number generalizes the Landau level number n: when the wave function does not 
depend on the second variable k, the iV-level classification reduces to usual Landau levels. 

In the case of anisotropic crystal, the rotational symmetry is broken. If, however, magnetic field 
is directed along a p-fold axis, one can see that iV-level classification reduces to classification in N 
over mod p. If magnetic field is directed in the basal plane of hexagonal crystal (along C2 axis) the 
above classification is reduced to the parity of N, as we have seen for the eigenfunctions (0) and 

In the absence of spin-orbital coupling spin projection on the field direction provides another 
quantum number, which distinguishes irreducible representations. If the spin-orbital interaction is 
strong, the mixing of singlet and triplet superconductivity is possible in a magnetic field. From the 
symmetry point of view, the space inversion P transforms basis functions of a particular Landau 
level through each other. Therefore, functions ( |64|) do not possess any parity under P. Hence at 
Hc2 the triplet states on which the same irreducible ray representation of the group Q is realized 
should arise simultaneously with the singlet states, although the pairing type occurred at H = 
has near Hc2 the dominant amplitude. 

In conclusion, one can apply the above classification scheme for all systems described by the 
symmetry group (0). Two-dimensional electron systems under a magnetic field are the most 
interesting among others. They were considered, e.g., in the works [■^^jTj^ where particular level 
operators analogous to Casimir operator (|62D were introduced. It is general symmetry properties 
of the system which make it possible such procedure in all cases: for a single electron, interacting 
electrons, electrons with spin-orbital coupling, or electron-phonon interaction. 

3.3 Kink in the Upper Critical Field 

The peculiar H-T phase diagram of superconducting UPt^ with two jumps in specific heat at 
temperatures and Tc2 and with the intersection of two critical fields Hi{T) and H2{T) at the 
kink point (see Fig. 1) is the main argument in favor of the anisotropic pairing in heavy fermions and 
of the multicomponent superconducting order parameter in UPts in particular p^]. The splitting 
of phase transition is believed to be explained by the closeness of critical temperatures of two 
superconducting states. The critical fields Hi(T) and H2(T) are interpreted as upper critical fields 
of these phases. 
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The explanation of the phase diagram of UPts has two different aspects: (i) what is the nature 
of the sphtting of the phase transition, and (ii) under what conditions the double superconducting 
phase transition leads to the phase diagram similar to that in Fig. 1. 

A comprehensive discussion of the first point can be found in [27,43,611 _ only enumerate the 
most popular models of Tc splitting in UPts. 

(1) Antiferromagnetic (AFM) symmetry breaking field models, (a) In some of these models p4,28,45,6ij 
it is assumed that the unconventional superconductivity in UPta is described by the two-component 
order parameter r] = {r]i,ri2) in accordance with (0) which corresponds to a 2D irreducible represen- 
tation of the point crystal group Dq, either Ei or E2. A weak AFM order with staggered moments 
in the basal plane lowers the symmetry group of the system from to D2 and therefore splits Tc. 
(b) A more complicated variant [^^] assumes a one-dimensional irreducible representation of Dq for 
superconducting order parameter but with additional degrees of freedom due to the triplet pairing 
and weak spin-orbit coupling. The splitting is due to the interaction of AFM moments with the 
total spin of the Cooper pairs. 

(2) The model of nearly isotropic d-wave pairing. In this picture f^] UPta is a d-wave isotropic su- 
perconductor with small splitting of critical temperature according to the irreducible representations 
of D^h group. Two of these representations Aig and Eig are responsible for the phase transitions in 

UPtg. 

(3) Accidental degeneracy models [^^'^S] unlike the other models do not explain the physical reasons 
for Tc splitting. They only suggest that the closeness of critical temperatures of two irreducible 
representations of Dq^ group is accidental. 

The answer on the second question can be given using the quantum numbers technique of Sec. 
3.2. The upper critical field corresponds to the maximal eigenfield of the linearized equations. The 
crossing of two eigenfields Hi(T) and H2(T) leads to the kink in Hc2(T) = ma.x{Hi(T) , H2(T)} . 
Using the general quantum mechanical rule that only symmetry different eigenlevels of the linear 
problem can cross, we came to the following conclusion: 

The kink is possible only if the eig enf unctions related to Hi{T) and H2{T) have different 
quantum numbers N or a. Otherwise the kink should he absent or at least smeared. 

Application of this rule to the cited above models of Tc splitting allows to verify easily their relevance 
to the phase diagram of UPts. In fact several models [i9'53,6i,83j satisfactory explain stability of the 
kink in Hc2{T) for different orientations of H, and a new experimental information is needed to 
select one of them. 

3.4 Anisotropy of the Upper Critical Field 

In conventional anisotropic superconductors the upper critical field can be calculated by taking into 
account the anisotropy of gradient terms in the GL functional. One should replace the gradient term 
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in (^) by KijD*^>* Dj^ (in the BCS theory, Kij is proportional to the inverse tensor of the effective 
masses of conduction electrons). The second rank tensor Kij has three independent components 
{K^x, Kyy, Kzz) for biaxial crystals, two components (i^n = K^z, K± = K^x = ^yy) for uniaxial 
crystals and only one component {K^x = Kyy = K^z) for cubic crystals. Therefore, in uniaxial 
crystals (e.g., with hexagonal or tetragonal structure) the angular dependence of Hc2{'&) has the 
form of an axisymmetric ellipsoid with the main axes: h}.2 = hca/2\e\K±, = hca/2\e\{K±K\iY^'^ 
(see Fig. 5). For a cubic crystal the upper critical field in the GL region near Tc is isotropic. 

Beyond the GL region the higher gradient terms slightly break the ellipsoidal form of the angular 
dependence of Hc2 and restore the anisotropy of Hc2 which coincides with the full point symmetry 
group of the crystal. However, this contribution to Hc2 is of the order of (1 — T/Tc)^, and it does not 
change the slope of the upper critical field = {dHc2/dT)\T=Tc- Recently different types of upper 
critical field anisotropy in the higher orders of (1 — T/Tc) have been discussed for conventional and 
unconventional superconductors in 

Turning back to the unconventional superconductivity, one should expect that anisotropy prop- 
erties of different from those for the s-wave superconductor due to a more complicated 
structure of the GL equations. The problem of anisotropy for different multicomponent order 
parameters for the crystal structures relevant to the heavy fermions was investigated in pi-H-'^^j^ 
The following unusual for s-wave superconductors features were found: 

• For the cubic crystals like UBeis (point group Oh), H'^2 should always have the cubic anisotropy 
for the different field orientations. 

• For the tetragonal crystals like CeCu2Si2, URu2Si2 (D^h), the square-like anisotropy of H'^2 is 
expected when the magnetic field lies in the basal plane (see Fig. 6) 

• For the hexagonal crystals like UPta, UPd2Al3, UNi2Al3 (D^h), there are no additional features 
of H'^2 anisotropy with respect to s-wave superconductors. H'^2 is isotropic in the basal plane 
and has a uniaxial outplane anisotropy. 

Note that these types of "unusual" anisotropy are connected to the symmetry of corresponding 
GL functionals. However, such anisotropy was found for neither UBei3 p], nor for CeCu2Si2 [^^]. 
Moreover, the upper critical field in CeCu2Si2 shows unexpected isotropy for all field directions in 
the linear regime near Tc [^^]. This discrepancy can be explained by the weak sensitivity of the 
electron-electron pairing interaction to the crystal anisotropy (e.g., it is the case for the p-wave 
pairing p^]). 

Another characteristic feature of the anisotropy of the upper critical field in the multicomponent 
superconductors is the nonelliptical, nonmonotonous behavior of Hc2{'&) for uniaxial crystals when 
the applied field changes its direction from the principal axis to the basal plane. This behavior is 
connected with the existence of several solutions with close critical fields for arbitrary directions of 
H rather than with the symmetry of GL equations p°]. The investigation of the GL equations for 
two-component model ([T5|) in [^'^°] shows that the angular dependence of H'^2 can have pronounced 
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minimum in a particular region of the phenomenological parameters C and K (see Fig. 7). Such 
behavior of Hc2{'d) was also found beyond GL approximation in [^'^], what confirms the universality 
of this feature. 

Unfortunately, the available experimental data on the anisotropy of the upper critical field in 
UPt3 [®^] are not clear enough to verify such a possibility. The observation of the nonmonotonous 
behavior of Hc2{'&) at uniaxial crystals CeCu2Si2, UPts, URu2Si2, UPd2Al3, or UNi2Al3 would be a 
good argument in favor of unusual superconductivity in heavy fermions. 

The striking feature of the superconductivity in UPta is the crossover in the anisotropy ratio 
H^/h12 with temperature, which is less than one near and greater than one at T = 0. Possible 
explanation of this experimental fact in terms of unconventional superconductivity is suggested in 

n 

At the end of this Subsection, we shall consider the quantum number classification for solution of 
the linearized anisotropic GL equations. When magnetic field is directed at an arbitrary angle with 
respect to the principal axis of a uniaxial system, the point symmetry group Dooh{Cooh) reduces to 
G' = {E,P,RU2x}, where x axis lies in the basal plane perpendicularly to H. At the first sight, 
there are no quantum numbers for solutions of linearized equations in this case because rotational 
symmetry about H is completely broken. 

However, the GL functional constructed for a single irreducible representation with a particular 
parity possesses an additional hidden symmetry. The action of space inversion P on a solution 
transitionally invariant along H reduces to the rotation by angle vr in the plane perpendicular to 
the magnetic field. Therefore, these solutions can be described by the parity of the Landau level 
number N, as in the case of H directed along the two-fold axis. For this reason, the vortex lattice 
in conventional superconductors has the form of the distorted triangular lattice with a two-fold 
symmetry at all intermediate directions of H. 

3.5 Modulation of Order Parameter Along Magnetic Field 

The appearance of superconducting phases modulated along field direction near Hc2 was discussed 
from the beginning of investigations of upper critical field in nontrivial superconductors [^^]. The 
discussion has been renewed last time in connection with the phase diagram of UPts [^'^'^^]. Here 
we formulate analytical approach to this problem on the example of two-component model. 

If the order parameter is allowed to vary along H, an additional quantum number, momentum 
Ph, appears. The reflection dh transforms functions with ph and —pn into each other, hence, they 
have the same eigenvalue. In Sec. 3.1 we have seen that for H ± c Eq. ( pOf ) could not be solved 
analytically for arbitrary px 7^ 0. Therefore, in this case we must verify the possibility that solutions 
with nonzero Px have smaller eigenvalues than 

It is obvious that solutions with px ^ Ijlu have larger eigenvalues because magnetic field may 
be considered as a perturbation to positively definite gradient terms in this case. We can 
check analytically also an opposite case oi Px <^ ^/Ih when instead nonzero Px can be treated as a 
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or 
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perturbation to eigenvalues (0) and (|52D. 

Let us consider Hamiltonian ( |50| ) as a sum of three parts HQ+H1 + H2, where Hq does not depend 
on px, Hi depends hnearly and H2 depends quadratically on p^. Note that Hi is an off-diagonal 
operator, therefore its contribution is of the order of pi and negative for the lowest eigenvalue. 
Eigenvalues and eigenfunctions of Hq are known: 

Ln = ^^^v^(2n + l), r/ = (/„,0), (65) 

nc 



Mn = ^''''™V ir(l + 2C) (2n + l) , rj={0,gn), (66) 

where fn and gn are introduced in (^TJ) and (^) wave functions of scaled Landau levels. They are 
characterized by the parity under reflection a^. and by the parity under rotation by angle vr about x. 
The variation of solutions along x breaks the classification according to the parity under reflection, 
but the parity under rotation by angle vr is still conserved. Then if we take rjx as a superposition of 
even Landau levels, riy must be a sum of wave functions from odd Landau levels and vice versa. 

In accordance with this symmetry classification, the second order perturbation theory gives the 
following corrections to the eigenvalues Lq and Mq: 

Loip.) = L, + Ki2,pI-KI,pI E (67) 



n=2fc+l ~ ^0 



Mo(;..) = M, + Kipl-Kl,pl E 
Substituting L„ and M„ to (^) and (|B5D, we find: 



2 



n=2fc+l ^" ~ ^0 



(68) 



Lo{p.)-Lo = Kiil + C)pl{l + e-2e\^il + qr±q''J^x 



1 



.2(1 + e) ((1 + e)i4k + 3) - VT^){1 + e + VT^] 



and the same expression for Mo{px) — Mq with the only change e —e. If C > 0, then Lq is 
the smallest eigenvalue for solutions uniform in the field direction. A lower estimate on Lq{px) is 
obtained neglecting the second term in square brackets. Then, using Y.T=o Jkal^^ = A — ; 
come to: 



Lo(p..) - Lo > ir,(l + Op^^ (1 + . - • (70) 

It may be easily checked that the right hand side of Eq. (|70|) is always positive for any C > 
(0 < £ < 1). The same is true for function M{p^) in the opposite case when —0.5 < C < 
(— 1 < 6 < 0). Thus, we may conclude that modulation of superconducting phase near Hc2 along 
H are unfavorable in both cases px <^ 1/Ih and px ^ ^/Ih- It is very improbable that such 
modulations will appear for intermediate values of Px- 
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While a solution with the lowest eigenvalue, e.g. Lq for C > 0, is stable against small nonhomo- 
geneous perturbations with ^ 0, this is not always true for the next level Mo(pxO- It follows from 
( |5BD that when Mq Li (i.e., C — 4) the corresponding denominator tends to zero, giving a large 
negative contribution to the coefficient at p^. Thus, at C = 4 the level Mq is certainly unstable 
against modulations with particular wavelength. The numerical investigation of Mq{px) shows that 
coefficient at p"^ changes sign from + to — at C = 0.946. 

A well developed modulated structure at C > 0.946 can not be analyzed on the basis of the 
second order perturbation theory. An infinite system of secular equations must be solved in this case. 
Variational calculations show that for C = 4 (e = 4/5) when Mq = 3 (in of units 2\e\HKi^/K /he) 
the next level after Lq = 1 has the eigenvalue Mq{p*) = 2.24 with the modulation wavelength 
A = 2-71 /pI ~ 4.7 Ih- The discussed instability of the second solution of the linearized GL equations 
is important for consideration of phase transitions inside mixed state (see Sec. 5). 

4 Vortex Lattices Near Hc2 

The next necessary step in the investigation of the superconducting phase transition at the upper 
critical field is determination of the vortex lattice form near Hc2- We describe in this Section 
nontrivial features of this procedure applied to the multi-component solutions of linear equations 
like SK-phase (||). 

The symmetry approach to second order phase transitions p^] allows us to make a few general 
conclusions before fulfilling calculations. First of all, the free energy near Hc2 has the form of the 
Landau expansion in powers of A. The structure of the order parameter A(k, r) is determined by the 
solution of the linear equations (|5^) and corresponds to the particular irreducible ray representation 
{N,a). The superconducting amplitude A is proportional to (iJc2 — HY^'^. The degeneracy of the 
order parameter is lifted by the forth order terms of the Landau functional. While the second order 
term is unique and can be always presented as |A^p, the forth order terms, generally, become 
nonlocal at low temperatures. For this reason, we consider only a restricted problem — investigation 
of nonlinear equations in GL approximation. 

We derive below the energy parameter for multicomponent GL functionals which generalies 
corresponding Abrikosov formula (Sec. 4.1). In Sec. 4.2 we discuss general procedure of finding 
the most stable lattice form and consider symmetry properties of the energy parameter under 
different choices of the lattice parametrasation. Sec. 4.3 is devoted to the search of the energetically 
favorable forms of the Abrikosov lattice for the SK-phase (^), which should appear near Hc2 in the 
two-component superconductor. Symmetry group of the vortex lattice is constructed in Sec. 4.4. 

4.1 Average Energy Density 

In this Subsection we construct the energy parameter for an unconventional superconductor, which 
determines the form of the vortex lattice. For this purpose we must take the order parameter in 
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an unchanged form as a solution of the hnearized GL equations with the largest critical field. The 
relative accuracy of this approximation is (1 — H/Hc2)- The vector potential Aq (rot Aq = Hc2), 
which enters the linearized GL equations, must be changed in order to include the diamagnetic field 
of superconducting currents: h = rot A = H + h^, rot = 47r/c ■ j^. Neglecting by the boundary 
effects, we come to the following functional: 

F=!^-^{\k.\} + {\A\'}-!f-!^, (71) 

An OTT 87r 

where (. . .) denotes spatial average. The forth order terms in (^) are written in a symbolical form, 
which implies the sum over all the possible invariants constructed from the components 

Thus the calculation of F reduces to the derivation of the diamagnetic field hg with the help of 
the Maxwell equations [^^] . For one-component superconductors, in which spatial dependence of the 
order parameter is restricted to zero Landau level functions, this field may be found analytically: 
hs = 8TTeK\^{r)\'^/nc. 

Since the calculation of hg for the multicomponent SK-phase becomes more difficult, we restrict 
ourselves to the case of k ^ 1. This allows us to neglect by the last term in (fflD, which is always 
of the order of (|A|^)//t^. Thus we should find only {\hs\). In the chosen approximation hg is the 
quadratic form of Ajv(k, r). The magnetic moment of the unit volume is a quantity invariant under 
action of the symmetry group (^31). Therefore, the average {\hs\) is the second order invariant 
constructed from A(k, r), and, consequently, {\hs\) oc (|AAr(k, r)p). The numerical factor in this 
relation does not depend on the particular structure of solution, because formally at if = the first 
term in (^) must be equal to — a(|Ap). Thus for large k the multicomponent GL energy is: 



Minimizing (^) by the amplitude of A, we get 

F = --(l 

where 

(|A|^) 



(74) 



(|A|2)2 

is the energetic parameter of GL model which generalies usual Abrikosov parameter for ordinary 
superconductors. The absolute minimum of 6 determines the geometrical form of the vortex lattice 
near Hc2- For conventional superconductors it is 

in accordance with the Abrikosov formula [^]. Vortex lattice form for axial phase (|48|) of the two- 
component model is determined by the same parameter (|75| ) with (3 changed on Substituting 
in (173) we obtain that for the SK-phase the energy parameter has the following form: 

(A(|/or + ^'1/21') + 2u;\(3^ + 2(32)\mf2\') 



JSK 



(76) 
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We discuss energy minimums of these parameters in the Sec. 4.3. 



4.2 Properties of Lattice Solutions 

Before search for the favorable lattice form we construct in this Section general periodic solutions 
on different Landau levels which then should be substituted in ([7^)- ([76|) . We will choose particular 
parametrisation for two-dimensional vortex lattices and discuss its symmetry properties which are 
helpful in further numerical calculations. 

The vortex lattice is defined by the discrete group of translations through a and b: 

fa A(k, r) = expitifi) A(k, r) , fb A(k, r) = exp{tif2) A(k, r) , (77) 



with arbitrary phase factors (fi and (p2- It follows from (^) that magnetic translations for this 
subgroup commute. Thus, due to (|57D any periodical solution possesses an integer number of flux 
quanta per unit cell p''''^]. 

We consider only lattices with one quantum of flux per unit cell. All distances are measured 
in the units of the magnetic length Ifj = hc/2\e\H, and the vector potential is chosen in the 
Landau gauge A = {—Hy, 0,0). Using the explicit form of the magnetic operator in this gauge 
Ta = exp{—iayx) exp(aV), it is possible to construct the lattice solution on the nth Landau level 
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which satisfies (^) with a = (a,0), b = (6 cos a, 6 sin a), ipi = tc, ip2 = 7r(p + 1). The area of the 
unit cell is fixed by the condition of flux quantization: ab sin a = 27c. The lattice form depends 
on the complex parameter t = p + ia = b/a exp{ia). The parameters = exp(7ri/3) and = i 
correspond to lattices with hexagonal and square symmetries. The periodic function /o(r | r) has 
one zero per unit cell located at (0, 0). 
The functions (|78D can be written as: 

/o(r I r) = ^1 f - t] exp f-^) , /„(r | r) oc (a+)"/o(r | r) , (79) 



where z = x + iy, ^i{z\t) is the Jacobi theta-function, and the bar denotes complex conjuga- 
tion. The change of the vector potential affects only the exponential factor in the expression for 
/o(r I r), leaving the r-dependence through ■(9i(z/a | r) for all gauges. The theta-functions are the 
only quasiperiodical analytical functions. This leads to the uniqueness of /„(r | r) with given peri- 
ods a, b and phases ipi, ip2- Substituting (^) into (0), we obtain all possible one-quantum vortex 
lattices near Hc2- 

In order to find which type of lattice is realized, one should substitute A7v(k, r|r) into the 
nonlinear energy functional, calculate the free energy F{t) for each r, and select r which corresponds 
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to the absolute minimum of F{t). Before calculations we establish some symmetry properties of 

function -F(t) according to p^], which generahzes ideas of p^]. 

The choice of the basis for two-dimensional Brave lattice is not unique. The same Brave lattice 

may be described by different r. Instead of a and b any pair of their integer combinations ka. + Ih 

and ma-|- nb, such that kn — Im — ±1, can be taken as basis vectors. This transformation changes 

f k l\ 

parameter r to r' = {m + nr)/{k + lr). The matrices , which do not change the orientation 

\m n J 

of the basis form the modular group SL{2,Z). Together with substitution a —a, b — >■ b they 
represent a complete set of possible parameterizations of the original Brave lattice. 

In a physical problem, the symmetry of the periodical system can be lower than the symmetry 
of the corresponding Brave lattice. The state is specified not only by the lattice, but also by 
additional parameters describing the structure and orientation of the each lattice site. For such 
systems transformation in the space of its parameters that consists only of a different choice of basis 
vectors changes the initial state. Thus we are dealing with the presence or absence of the symmetry 
under the action of the group SL(2, Z) x Z2 in the complex plane of parameter r. 

Functions /„(r | r) (or AAr(k, r | r)) are uniquely determined by their lattice, and, consequently, 
they should be symmetrical under SL[2, Z) XZ2 group. This statement can be verified by establishing 
the following properties of Ajv(k, r | r) as a function of r (see proof in p^]): 

A^(k,r|T + l) = A^(k,r|T) , (80) 
AAr(k, r I - l/r) = e''^'^ exp 

where exp{iip) = ^t/t*, and 

A^(k,r| -T*) = i?[/2^ A^(k,r|T) . (82) 

Substituting Aiv(k, r | r) to an arbitrary functional F invariant under group Q, we obtain energy 
F(t), which is not changed under SL{2,Z) x Z2 transformations: 




The fundamental region of SL{2, Z) x Z2 group in complex plane contains all different one-quantum 
lattices. A specific example of the fundamental region is shown in Fig. 8. Their vertices Th and 
Ts are always stationary points of F{t). While at the point Tg functional can achieve all types of 
extremums, in Th only a local maximum or a local minimum are possible. The absolute minimum 
of F{t) can be located, of course, not only at these symmetrical points, but somewhere in the 
fundamental region. Moreover, such exotic behavior is characteristic for unconventional pairing 
(see Sec. 4.3). 

It is interesting that similar theoretical constructions have appeared independently in the theory 
of Weinberg-Salam model for the electroweak interactions [^^'^^]. Vacuum reconstruction in high 




Ajv(k,r I r) , 



(81) 
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magnetic fields in this theory is described by energy functional on complex wave function \E' with 
nonlocal forth order term arising from coupling with additional bosonic field. This nonlocality in 
the forth order terms could results in its turn in exotic forms of lattice solutions. 

4.3 Vortex Lattice in Two- Component Model 

The explicit dependence of the energy parameter for one-component superconductors from r can 
be obtained after substitution of /o(r | r) into 

5o(p, a) = py^ cos(27rpmn) e--(™'+"') . (84) 

m,n 

Here summation is fulfilled over all integer m and n. The absolute minimum of the function 6o{p, o") 
is achieved at the point of the fundamental region which corresponds to the regular triangular 
lattice. 

The energy parameter of the SK-phase (^) can be presented as: 

^sk{p, ct) = -^^^^^ (5oo(p, cx) + 2uj\l + 2/3)5o2(p, a) + uj%2{p, a)) , (85) 

where parameter u is given by (^) and P = ^2/^1- Explicit expressions for Spg{p, cr) can be found in 
Appendix. Positions of the absolute minimums are different for each function 6pq{p,a). Therefore, 
variation of phenomenological parameters (3 and u leads to the rich phase diagram for the vortex 
lattice form near Hc2- Note that, according to (^91), the following restriction < cu < I/V2 is 
fulfilled in the region of the existence of SK-phase near Hc2- 

The numerical investigation of 6sk{p,(x) was done first in [^^]. But the phase diagram found 
there contradicts in fact general requirements of the phenomenological theory, because it predicts 
two coinciding phase transitions lines between triangular and square lattices and between square 
and rectangular lattices. The correct phase diagram for the position of the absolute minimum of 
function 6sK{p,cr) versus phenomenological parameters (3 and 00"^ is presented in Fig. 9. There are 
different regions where the energetically favorable solution corresponds to the regular triangular, 
square, rectangular, and distorted triangular lattices. Phase transitions lines in Fig. 9 between 
square and the rectangular lattices and between square and distorted triangular lattice IV are of 
the second order, all other transitions are of the first order. For negative (3 the only stable form of 
the vortex lattice is a triangular one. 

The characteristic scale of the difference between various extremums of Ssk{p, o") is of the order 
of 10~'^-10~^ compared to 10"^ for 6o{p,a). Note that the point (/5,u;^) = (0.5,0.1) predicted by 
BCS theory for Eig representation is situated in the region of stability of the hexagonal lattice. It is 
interesting that for Ssk{p, cr) the absolute minimum always lies on the boundary of the fundamental 
region. 

In Sec. 4.2 we have argued that the one-quantum vortex lattice with hexagonal symmetry always 
corresponds to the extremum of the energy functional for arbitrary superconducting order parameter 



75). This calculation was done in 
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A7v(k, r). Now we can see from Fig. 9 that corresponding stability region of the hexagonal lattice 
for the SK-phase is rather large. 



4.4 Symmetry Group of the Vortex Lattice 

To complete the discussion of vortex lattices near Hc2 we should write explicitly corresponding 
symmetry groups, which include magnetic translations and rotations ( ^5l) in combinations with 
nontrivial phase factors. 

The transformation rules of Ajv(k, r) near ifc2 under magnetic translations by the lattice periods 
are defined by (^) and ([781). Vortex lattices corresponding to different phases ipi and go over 



into each other under magnetic translations which do not coincide with any integer combination 
of the basis vectors. Displaced solutions TpA7v(k, r | r) are transformed like Bloch wave functions 
with wave vectors lying in the Brillouin zone of the background lattice A7v(k, r | r) \^'^]. 

We are interested mainly in the transformation properties of the triangular lattice, since it has 
a particular region of stability for arbitrary N. Using identity — 1/r/i = Th one can obtain from 
(p0|) -(^2l) following properties under rotation by angle 7r/3 about arbitrary vortex axis F: 

A;v,r(k, v\Th) = e-^-(i-^)/3^;v,r(k, r | , (86) 

and under time-reversal together with rotation by angle vr about an axis perpendicular to the 
direction of magnetic field: 

RU2y AN,r (k, r | n) = AN,r (k, r | r/,) . (87) 

Combining the elements listed above, we can write the symmetry group of the vortex lattice as 

6^r = {e*"Ta, e'^/^n, e'^^'-^y^L^/,, aa^, RU2y} . (88) 

All other elements are formed as proper combinations of this basis. For example, the supercon- 
ducting order parameter AArj(k, r | r/^) also possesses symmetries under rotations by angle 27r/3 
about an axis passing through the center of the triangle K and on angle vr about the middle of the 
triangle side M (see Fig. 10). Considering two other displaced vortex lattices Ajv ;^(k, r | t/j) and 
A7v,M(k, r I Th), whose points K and M coincide with the center of coordinate F, we can establish 
their symmetries with respect to F: 

Qk = {e-'^'/'f^, e-*-/6fb, e-2*-^/3L2V3, (rat, e'^RU2y}. (89) 

QM = {e'^f^, e-^'^/^fb, e^'^^L^, aan, e'^RU2y]. (90) 

The symmetry properties of square lattice A7v r(k, r | r^) and of all others from Sec. 4.3 can be 
investigated by the similar way. 

We have derived transformation properties (p8D-(pO|) using explicit form of A7v(k, r | r) which is 
valid only for a constant magnetic field inside a sample, i.e., when H Hc2- Below Hc2 the magnetic 
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field is modulated by superconducting currents. However, Eqs. ( pSj) allows introduction of discrete 
magnetic translations and rotations in this case too. Therefore, the mixed state is invariant under 
symmetry group ( |55D in the finite region H* < H < Hc2 until a new phase transition occurs. 

5 Phase Transitions in the Mixed State 

Perhaps the most unusual feature of the mixed state of an unconventional superconductor is the 
possibility of the phase transitions between different vortex lattices. 

As we have mentioned in Sec. 1.2, in usual superconductors the character of vortex interaction 
strongly depends on the distance between vortices. In the vicinity of the lower critical field Hd 
vortex lattice is formed through the interactions between nearest neighbors. At fields H > He 
vortex interacts with all neighbors within circle of the radius A. If, however, H <^ Hc2 the average 
distance between vortices is much greater than coherence length ^ and the intersection of cores 
of different vortices is negligible. In this case only long range hydrodynamics repulsion participates 
in the lattice formation. 

As vortices in unconventional superconductors have the same structure at distances r > C,, all 
above arguments are valid for them too. The cores intersection becomes important for intervortex 
distance ~ max{^, ^}. Therefore all deviations from the usual Abrikosov theory of vortex lattice 
formation can occur only in the field region near the upper critical field. The convenient tool for 
investigation of vortex lattices in that region is symmetry approach developed in Sees. 3 and 4. 

In Sec. 5.1 we formulate analytical approach to the investigation of phase transitions near Hc2 
and show that the closeness of two critical (eigen) fields in the linear problem always leads to 
additional phase transitions inside mixed state. Sec. 5.2 is devoted to the symmetry consideration 
of structural phase transitions in the vortex lattices. Finally, in Sec. 5.3 we briefly discuss possible 
second order phase transitions below upper critical field on the examples of two-component GL 
functional (H) and accidental degeneracy models. 

5.1 General Approach 

The space distribution of the order parameter A(k, r) at a given temperature T and an external mag- 
netic field H {Hci < H < Hc2) is defined by the minimum of the energy functional F{A(k, r), H, T}. 
To minimize F, one needs to solve the nonlinear (GL) equations: SF/SA* = 0, 6F/6A = which 
can be symbolically written as: 

L{A(k,r),/7,T} = 0. (91) 

Their solution near Hc2{T) is the Abrikosov lattice with a particular symmetry group Qr- 

Suppose that at some critical field H*(T) the second order structural phase transition in the 
vortex lattice does take place. Speaking in terms of the phase transition theory, H*{T) is the 
bifurcation line below which Eq. ( pi]) has two different solutions: the old one Ai(k, r) invariant 



34 



under the group Qr and the new one Ai(k, r) + A2(k, r) with a lower symmetry Q. At H < H*{T) 
the energy functional changes the character of the extremum for solution Ai(k, r) from the absolute 
minimum to the local maximum, whereas Ai(k, r) + A2(k, r) corresponds to its minimum. The 
small perturbation A2(k, r) belongs to some nonunit irreducible representation of the group Qy- 
Since the phase transition is of the second order, the amplitude of A2(k, r) at if < H* is small as 
[H* — HY^"^. Linearizing (^) in the small value A2, one obtains the equation for the H*(T) as a 
condition on the magnetic field at which the linear eigenvalue problem 

£{Ai(k,r),ii*,r}A2(k,r) = (92) 

has a solution. The linear operator £. depends on the space distribution of the order parameter 
Ai(k, r) which in its turn is the solution of (|9T|) taken at the phase transition line H*{T). 

To specify the influence of the background solution Ai(k, r) on the appearance of A2(k, r), we 
separate the part Cq from the operator C which does not depend on A and present ( P^ in the 
following form: 

Co{H*iT)} A2(k, r) + £i{Ai(k, r), H*{T)} A^iK r) = . (93) 

In the vicinity of ifc2 the second term in (p3| ) is small and can be considered as a perturbation to 
the equation 

Co{H,T}A{k,T)=0, (94) 

which in fact coincides with (^41). The last equation has different critical fields Hi{T), and the 
maximal of them Hi{T) corresponds to Hc2{T). Let us assume that there exists another eigenfield 
H2{T) of Eq. ( p^ ) which is sufficiently close to Hi{T). Then one can derive similarly to ( |72D the 
following two-order parameter functional: 

F= - a,(l-^) (|Ai|^) + (lAil^) -a,(l- -|) {\A,f) + {\A,\') + {AjAj) , (95) 

which represents free energy for fields close to ifc2- 

Due to the infinite degeneracy of both order parameters, finding all extremums of the energy 
functional (^) is a rather complicated problem even in the case when spatial dependencies of both 
Ai and A2 are described by zero Landau level functions. However, without interaction between Ai 
and A2 the problem of the phase transitions described by functional ( p5|) becomes trivial. There 
will be two successive second order phase transitions at H = Hi and H = H2. Turning on of the 
interaction term (A^A2) affects only the lower transition. At least in the finite range of parameters 
where interaction term is small, the type of the lower transition is not changed and the critical field 
H2{T) is slightly renormalized to H*{T), which corresponds to the structural phase transition in 
the vortex lattice. This case is considered in the two following Subsections. 

If other eigenfields of (0) are far from Hc2{T), this analytic approach is not applicable, because 
Eq. (p3|) can not be reduced to (0) due to the significant admixture of all other Landau levels. 
However, the continuous change of parameters of the system, which move H2 from Hi to lower 
fields, raises a hope that in some range of parameters the phase transition in the mixed state still 
exists. 
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5.2 Structural Phase Transitions 



According to the Landau theory, if the second order phase transition in the vortex lattice occurs, 
the initial symmetry cannot be broken to an arbitrary subgroup of Qr- The residual symmetry 
group ^ C is defined as the symmetry of the new lattice Ai(k, r) + A2(k, r). To classify the 
possible ways of symmetry breaking in the vortex lattice Ai(k, r) one should (i) find all irreducible 
representations of the group Qr, (ii) decompose the new order parameter A2(k, r) over the basis 
functions of a particular representation, (iii) write the corresponding Landau functional for chosen 
representation, and (iv) determine all its minimums with their residual symmetry groups Q. 

We do not intend to develop here this general approach. Staying on the point of view of the 
above Subsection, we only show how this symmetry analysis must be applied in the particular case 
of the energy functional (pSl). 

Previously we have shown that the appearance of the phase transition in the vortex lattice at 
H = H*(T) should be connected to some eigenlevel of (0). The quantum numbers (A^2,o"2) of 
A2(k, r) are different from quantum numbers (A^i,cri) of Ai(k, r). Therefore, these solutions are 
orthogonal: {AIA2) = 0. Moreover, the existence of a structural phase transition at H = H*{T) on 
a nonunit irreducible representation of the group Qr also means that (|AipA]^A2) = 0. Otherwise 
nuclei of A2 would arise at if = Hc2{T) simultaneously with the vortex lattice Ai, as is the case 
for all Landau level functions with n (mod 6) = in conventional superconductors p^] . 

Irreducible representations of two-dimensional space group Qr are characterized by the wave 
vector q lying in the Brillouin zone of the vortex lattice. To obtain a set of corresponding Bloch 
functions from those belonging to the Landau level N2, we should construct one-quantum lattices 
A2 of the same form as Ai. Then, as we have discussed in Sec. 4.4, there exists one-to-one corre- 
spondence between the displacement of periodical solutions A2 in the unit cell of Ai and the wave 
vector q in the Brillouin zone. 

According to the Lifshitz criteria, only symmetrical points of the Brillouin zone are impor- 
tant for the structural phase transitions in a crystal lattice [^^]. For the regular triangular lattice 
these are the center of the hexagon, its vertex, and the center of the hexagon side. These points 
of the Brillouin zone correspond to the following displacements: pp = 0; Px k' — ^^(a + b)/3; 
Pm,m',m" = a/2,b/2, (a — b)/2. The existence of equivalent displacements is a consequence of 
rotational symmetry around F. In the case of a distorted triangular lattice, e.g. when H ± c, this 
degeneracy is lifted. 

The shape of the stable lattice structure Ai + A2 on the transition line H*{T) is determined by 
interaction in the quartic terms: (|Aip|A2p) and, if it exists, (A*^A2 + c. c. ). Depending on the 
sign of these interactions, all relative locations between Ai and A2: F-F, T-K, F-M are possible. 

The residual symmetry group of the mixed state ai H < H* is determined by 

g = g^^n g^' (96) 

with appropriate Q = T, K or M. 



36 



If the vortices of the new lattice coincide with those of Ai(k, r), only rotational symmetry is 
broken below H*. This implies the distortion of the initial hexagonal lattice. 

If new vortices appear between the old ones in M position, the operator e~*'^^Tb transforms 
A^'"^ + Ag*^'' into Aj^'"'' — h.Q^^'^ (see (|^) and (0))- In this case the unit cell of the whole lattice is 
two times larger than separate cells of lattices Ai and A2. It follows from the same arguments that 
+ Af ) carries three quanta of the magnetic flux per unit cell. Such period multiplications can 
be either complete or partial when elementary translations Ta or Tb belong to Q in combination 
with proper rotations or reflections. 



5.3 Examples 

5.3.1 Closeness of two eigenfields 

For s-wave superconductors the upper critical field always corresponds to the zero Landau level. 
The first Landau level has a critical field three times as small. For multicomponent superconductor 
it is possible situation when eigenfields of different phases lie sufficiently close to each other. This 
leads to structural phase transitions inside the mixed state. 

To illustrate this case, we use the example of the two-component superconducting order param- 
eter. As we have seen in Sec. 3.1, depending on the constants C and D two different phases can 
arise near Hc2 for H || c. For example, if \sk < ^a, the SK-phase with the generalized Landau level 
number Ni = +1 appears at Hc2 = Hsk- If, however, Xa is slightly larger than Xsk, the critical 
field Ha for the axial phase with Ai"2 = — 1 is close to Hc2- When the regular triangular lattice is 
favorable for SK-phase, the nuclei of axial phase do not appear simultaneously with SK-phase at 
Hc2 because triangular lattices for phases with different are invariant under different hexagonal 
groups (il)- 

To investigate possible phase transitions in the vicinity of -ffc2 one should use energy functional 
(^51) with ai = a2i nonlinear terms (^) and (|76D , and the interaction term of the form (k ^ 1): 

(A^Ak^) = 2Ac.2(i/oiv2"^r) + (A+2/52)(i/o"n/o^''r) • (97) 

According to the results of p^], the instability aX H = H* (minimum of (^)) corresponds to the 
appearance of T-K structure if (32/ Pi > —0.5 — 0.12^7^ with the residual symmetry group 

^ = {e*-/^fa+b, e-^'^/^Taa-b, e'^^'^/^^L^^g, aa,, RU2y} . (98) 

If —0.5 — 0.12u;^ > P2/ Pi > —0.5 — 0.35 cj^ new vortices are located at M positions, and 

g = {e'^f^, e'-fab, aah, RU2y} . (99) 

If —0.5 — 0.41 a;^ > P2/P1 > —1 two lattices Ask and are not displaced. Their symmetry group 
is 

^={e^'^ra, e^'/'Tb, L,, aa^, RU2y} . (100) 
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For parameters — O-SScj^ > /52//5i + 0.5 > —QAl omega? the minimum is achieved for nonsymmetric 
displacements between M and F points. 

One can see that the period multiphcation in the first two cases is complete. That is, e.g., the 
diamagnetic field of superconducting currents hg is periodical only under translations from (^) and 



The phase transition at if = H* leads to the phase diagram of the two-component GL functional 
presented schematically in Fig. 4. 

For H ± c there are also two phases (|5l|) and (jS^) with close values of critical fields whenever 



C is small. Since these phases differ in the quantum number a and in the parity of A^, their mixing 
is prevented at iic2- Possible structural phase transitions in this case were partially considered in 
p6,82j r^YiQ complete analysis over the entire region of phenomenological parameters must take into 
account the possibility of the appearance of phases modulated along H (see Sec. 3.5), and has not 
been done yet. Note that the available experimental data on the vortex lattice restructuring at 
H = H*{T) in UPta correspond to this orientation of magnetic field relative to the crystal axes 

The example in this Subsection shows that for any pairing type (even in an anisotropic crystal 
with the magnetic field directed along the rotational axis) in some range of coefficients of energy 
functional on the H-T phase diagram there exist a phase transition line at Hci{T) < H*{T) < 
Hc2{T). How large is this region is the question of numerical calculations. 



5.3.2 Intersection of critical fields belonging to the different pairing types 

Let us assume that at H = two different superconducting states with close critical temperatures 
Ta and (T^ > T^) do occur. Phases a and b may belong to the different pairing types {s + p 
OT s + d) in an isotropic metal or to the two irreducible representations of the point group G in 
anisotropic crystal. Here we will not discuss the physical reasons for this accidental closeness of Ta 
and Tft but restrict ourselves to the question: what is the H-T phase diagram of the system in this 
case? 

The coexistence of two superconducting phases in magnetic field is described by some GL energy 
functional: 

F = Fa{Aa} + FbiAb} + Fab{Aa, ^b} , (101) 

which takes the form of (p5|) near the upper critical field. The solution of the problem of H-T 
phase diagram is trivial if there is no interaction Fab{^a, ^t}'- if will be the sum of two independent 
diagrams for the phases A^, Ab with the upper critical fields Ha, Hb, and under the condition 
dHb/dT > dHa/dT it will resemble the phase diagram of UPta (Fig. 11). 

Being described by different irreducible representations of the group G at H = 0, the phases A^, 
Ab can also belong to the eigensolutions of linearized equations of a different symmetry. Following 
Sec. 3.2, we should attribute the different quantum numbers {Na,<Ja), {Nb,<Jb) to the critical fields 
Ha, Hb- The interaction term Fab{Aa, Ab} does not change the symmetry classification of the 
eigensolutions of the linear problem. It leads to the renormalization of Ha{T) and Hb{T) curves 
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inside the mixed state but does not smear the kink in Hc2{T) = ma.x{Ha{T), Hfy{T)}. Phase 
diagrams for different combinations of irreducible representations for and based on Landau 
expension near Hc2 ( P^D have been partially considered in ps.ig.ssj^ Note, that without mixing of 
Aa and A;, in the higher order gradient terms the model ( |101| ) exhibit only partial breaking of 
translational symmetry below H*{T). 

The field intersection is also possible in the above example of 2D order parameter. The small 
uniaxial anisotropy in the basal plane {\i]x\'^ — 1%^) splits Tc into: for $y(k) phase and T~ for 
$3;(k) phase. If the magnetic field is directed along x, there exist an intersection of two critical 
fields: Hy{T) ~ (T+ - T) / {Ki2zKif/'^ and H'^iT) ~ {T~ - T)/{KiK^y/\ In fact it is the first 
model proposed for the phase diagram of UPts [^■*''^^] . 

Concluding this Section, we want to emphasize that the quantum number technique is very 
important for establishing the nature of the phase transition at H*{T) as well as for investigation 
of the stability of the kink. The exact sequence of the phase transitions for arbitrary strength of 
the interaction term in the energy functional (^) is still unknown. By numerical calculations in [^^] 
it was shown the possibility of a phase transition of the first order for large strength of interaction, 
which occurs before the discussed second order structural transitions. 

To summarize, we have reviewed the equilibrium properties of the mixed state of superconductors 
with the multicomponent order parameter. The main feature of unconventional pairing is the 
reduced symmetry of the superconducting state in comparison to the normal one. This results in a 
large variety of unusual properties such as power law temperature dependencies of thermodynamic 
characteristics, domain structure due to the presence of degenerate states, the angular dependence of 
the Josephson current, and others, which were considered in details in previous reviews [22'65j_ 
perhaps the most pronounced breakdown of various symmetries occurs when vortices penetrates 
inside the volume of a multicomponent superconductor. We have considered nonaxisymmetric 
vortices and flux lattices of different forms. Structural phase transitions in the Abrikosov lattice 
are the main qualitatively new feature of unconventional superconductivity which, in our opinion, 
should be primarily used for establishing unusual superconductivity experimentally. 
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APPENDIX 

The energy parameter (|T^) for unconventional superconductors is constructed in general case from 



the repeated blocks which are integrals of periodic functions from different Landau levels. In the 
case of the SK-phase they are 

^oo{p,(^) = MiX = E cos(27rpmn) e-"'^(™'+"') , (102) 



my 

{\J0\ ){\J2\ ) ^ m,n \4 

+ T,^a\m'-n'f\, (103) 



^22(P, a) = 1^ = ^ E cos(2vrpmn) e---(™^+"^) (11 - ^vra(m^ + n^) 

\\J2\ I ^ m,n \IK) Z 

+ 37rV2(m2 + n2)2 + ^^T,^a\m' - - (104) 
- 67rV3(n2 + m2)(n2 - + 7rV^(n2 - m^f^ . 

Behavior of these functions on different lines in the fundamental region is shown in Fig. 12. The 
corresponding absolute minimums are located at the following points: for (5oo(p, cr) at p = 0.5, 
(7 = for 5o2(p, 0") at p = 0, cr = 2.6, for 522(p, a) at p = 0.5, a = V3/2. 
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Fig. 3. Phase diagram of axisymmetry instability at the vortex core in the two-component model 
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Fig. 4. The eigen critical fields and schematic H-T phase diagram for H || c in the two-component 
model when D < C^/(l -|- C). For D > C^/ (1 -|- C) the order of eigenfields is reversed. 

Fig. 5. Outplane anisotropy of Hc2{'&) in uniaxial s-wave superconductors; K\\ : = (a) 1 : 2, 
(b) 2 : 1, (c) 4 : 1. 

Fig. 6. Anisotropy of Hc2{<f) in the basal plane for Ei irreducible representation of group, Kr, 
is the coefficient at the "tetragonal" invariant in gradient terms |Tl[] and K2 = K3 = Ki; Ki 
: i^s = (a) 2 : 1, (b) 1 : 1, (c) 1 : 2. 

Fig. 7. Angular dependence of Hc2{'&) between c-axis and the basal plane in the two-component 
model (5), e = C/(l + C),K = K^/Ki- (a) e = 0.3, K = 0.8, (b) e = 0.3, K = 1, (c) e = 0.3, 
K = 1.2, (d) e = 0.2, K = 1. 

Fig. 8. The fundamental region of the SL{2, Z) x Z2 group in the complex plane of the parameter 
T = p + ia. 

Fig. 9. Phase diagram of the vortex lattice form near Hc2 for SK-phase. Region I corresponds to the 
regular triangular lattice, region II — to the square lattice, region III — to the rectangular 
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a > 60° {a < 30°) respectively. Point corresponds to the weak-coupling parameters. Details 
of the phase boundaries intersections are shown on (b). 

Fig. 10. Symmetry elements of the triangular vortex lattice. 

Fig. 11. Schematic H-T phase diagram of the accidental degeneracy model. 

Fig. 12. Behavior of functions (a) 5oo(p, o"), (b) 6o2{p,(y), (c) 622{p,cr) on lines p = and p = 0.5 
in the fundamental region. 



45 



